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Abstract 

These notes form an introduction to cosmology with special emphasis on large scale structure, 
the cmb anisotropy and inflation. In some places a basic familiarity with particle physics is 
assumed, but otherwise no special knowledge is needed. Most of the material in the first two 
sections can be found in several texts, except that the discussion of dark matter and the cos- 
mological constant is more up to date. Most of that in the remaining sections can be found in a 
review of structure formation and inflation done with Andrew Liddle, which describes original 
work by various authors including ourselves and Ewan Stewart. The reader is referred to these 
works for more detail, and a very complete list of references. 
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1 The Recent Universe 



These notes form an introduction to cosmofogy with special emphasis on large scale structure, the 
cmb anisotropy and inflation. In some places a basic familiarity with particle physics is assumed, but 
otherwise no special knowledge is needed. Most of the material in the first two sections can be found 
in several texts Q , except that the discussion of dark matter and the cosmological constant is more 
up to date. Most of that in the remaining sections can be found in a review of structure formation and 
inflation done with Andrew Liddle which describes original work by various authors including 
ourselves and Ewan Stewart. The reader is referred to these works for more detail, and a very 
complete list of references. 

According to current thinking, the history of the observable universe broadly divides into three 
stages. First there is an inflationary era, when the energy density is dominated by the potential of 
a scalar field. Then there is a radiation dominated era when the energy density is dominated by 
relativistic particles, which are called 'radiation' by cosmologists. Finally, lasting till the present 
epoch, there is a matter dominated era when the energy density is dominated by the mass of non- 
relativistic particles, which are called 'matter'. This first lecture concerns the latter part of the 
matter dominated era, and in particular the present epoch. 

Unless the contrary is implied by the specified units, I set h = c = = ^■ 

1.1 How one decides what is true 

Before beginning let us address an important issue. More than most scientists, cosmologists tend 
to be asked 'how do you know?' The answer is that one proceeds in the same as in other areas of 
physics and related sciences. One starts with a model, which is a statement about the nature of 
physical reality, plus a theory formulated in terms of equations. Then one calculates ('predicts') the 
outcome of observations and experiments that have been performed, or better of ones that are about 
to be performed. If the agreement is impressive enough everyone agrees that the model is correct. 

This sequence cannot be reversed, that is one cannot rigorously work back to deduce the model 
from observation. To put it the other way round, a model can strictly speaking only be falsified, 
not verified. Logically, two completely different models could both agree with observation but that 
doesn't seem to happen, though of course an old model can be seen to be a limiting cases of a better 
one. 

So the answer to the question 'how do you know' begins with the admission that a guess has been 
made. At any one time, though, parts of the picture are universally accepted as correct, because 
they agree with such a huge body of observation. Other parts of the picture are just coming into 
focus, and agreement about the truth lies in the future. The scope of the agreed part is dictated 
not by logic, but by consensus. 

Except for the mathematical theory, all of this is not really special to science, but is rather 
just the way that human beings go about things. For instance, at any given epoch during the 
centuries that the geography of the earth's surface was being established in the West, parts of it 
were universally accepted because the accounts of many people agreed (or sometimes because just 
one of two surveyors had done what looked like a good job). At the same time, parts of it were still 
in dispute. The portion of the earth's surface about which there was agreement was dictated not by 
logic, but by consensus. 

Because our knowledge has been arrived at by comparing a large body of ideas with a large 
body of observation, it is not usually possible to point to a particular observation as the reason for 
believing a particular fact. For this reason, as well as because of the introductory nature of the 
lectures, I will often state a fact without giving specific justification. On the other hand, by the 
time we have finished it will hopefully be clear that the picture I put forward leads to agreement 
with an impressive body of observation, and that this picture could hardly be subjected to drastic 
alteration without spoiling this agreement. 

1.2 The geography of the observable universe 

Astronomers use the unit 1 pc= 3.26 light years= 3.09 x 10^^ metres. For cosmology the megaparsec, 
1 Mpc — 10^ pc is appropriate. The unit of mass is the solar mass Mq — 1.99 x 10'^^ g. 

Stars have mass in the range roughly 1 to IOMq. They are found only in galaxies may be regarded 
as the basic building blocks of the universe, with masses ranging from maybe lO^M© (dwarf galaxies) 
to IO^^Mq (large galaxies like our own). A galaxy typically has a luminous centre containing nearly 



all of the stars, and a dark halo of unknown composition which extends of order 10 times as far and 
contains of order 10 times as much mass. In many (all?) cases there is a nucleus consisting of a 
black hole which is gobbling up stars and gas and emitting radiation. If the nucleus is the dominant 
feature it is called an AGN (active galactic nucleus), examples being Seyfert galaxies and quasars. 
(Astronomers tend to using different names for different examples of the same type object. The 
reason, of course, is that they don't know at first that they are examples of the same type of object.) 

In round figures, large galaxies like our own have a size of .1 Mpc (including the dark halo) and are 
of order 1 Mpc apart. Many galaxies belong to gravitationally bound clusters containing from two to 
~ 1000 galaxies. Small clusters are usually called groups. Big clusters, of order 10 Mpc in size, are 
the biggest gravitationally bound objects in the universe. There do exists, though, 'superclusters' 
with size of order 100 Mpc, which are conglomerations of galaxies and clusters representing regions of 
space which have higher than average density. Presumably they will become gravitationally bound 
at some time in the future. On the scale of 100 Mpc there also seem to be sheetlike and filamentary 
structures, as well as voids containing hardly any galaxies. 

On scales bigger than 100 Mpc the distribution of matter in the universe is known to be very 
homogeneous, both from direct observation of the galaxies and from the isotropy of the microwave 
background. To be precise, if one throws down at random a sphere with radius R and measures 
its mass M, then the rms variation AM/M is a decreasing function of R, which is of order 1 at 
R = 10 Mpc and of order .1 at i? = 100 Mpc 

The biggest distance we can observe is of order 10^ Mpc, the distance that light has travelled 
since the Big Bang. The sphere around us with this radius is thus the observable universe. As far 
as we can tell, distant parts of the observable universe are much like our own part, with the same 
kind of galaxies, clusters and so on. From the fact that the microwave background anisotropy AT/T 
is of order 10~^, one can deduce that AM/M < 10~^ on scales R comparable with the size of the 
observable universe. 

What's beyond the observable universe? 

What is the universe like outside the patch we observe? Since the universe is very homogeneous 
and isotropic in the observed patch, it is reasonable to suppose that it remains so for at least a few 
orders of magnitude further. But what about happens after that? 

Until the last decade or so, the prevailing assumption seems to have been that the entire universe 
is homogeneous and isotropic, like the bit that we observe. To avoid the embarrassing question of 
what would then be beyond the edge of the 'entire' universe, this view requires that one believes in 
a 'closed' universe, which is allowed in the context of general relativity. Such a universe is the three 
dimensional analogue of a sphere, and as I discuss later on this lecture its spatial 'curvature' could 
be detectable observationally. 

More recently, the favoured view has been that as one examines larger and larger scales, the 
universe becomes steadily more inhomogeneous and anisotropic Q. Going the other way, if one 
takes an overall view of a patch of the universe many orders of magnitude bigger than our own, 
it has a non-fractal nature; as one goes down to smaller and smaller scales, it looks more and 
more homogeneities, as opposed to the 'fractal' situation where new inhomogeneities would reveal 
themselves at every step. On this view, the observable universe is extremely homogeneous because it 
is a small part of a much large patch which is roughly homogeneous and isotropic. As I shall discuss 
in Lecture 6, this viewpoint is normally discussed within the context of inflation which indeed makes 
it useful and rather compelling, but one should understand that it is a separate hypothesis, not 
related per se to the inflationary hypothesis. 

On this second view, there can be yet another hierarchy as one moves beyond the roughly homo- 
geneous and isotropic patch. One could encounter regions where the universe is radically different, 
for instance contracting instead of expanding, or with different laws of physics, corresponding to 
different solutions of some super-theory. One would presumably still like the 'entire' universe to 
be closed, but on this second view it would be analogous to the surface of an extremely deformed 
'sphere' so that no hint of the global behaviour would be evident from observations of our own 
region. 

1.3 The isotropically expanding universe 

On the large scales where it is homogeneous, the universe is also expanding isotropically. That is, the 
distance between any pair of galaxies separated by more than 100 Mpc is proportional to a universal 



scale factor a{t), the same for every pair. The Hubble parameter is defined hy H = a/a, where the 
dot is the time derivative. 

A subscript is generaUy used to denote quantities evaluated at the present epoch. It is con- 
venient to set ao = 1, so that a{t) is simply the size of any comoving region (one moving with the 
galaxies) relative to its present size. The present value of H, denoted by Hq is called the Hubble 
constant. It is traditionally measured by observing the redshift z = A A/ A of galaxies receding from 
us with velocity u <C 1. The velocity of such a galaxy is given hy v = Hr, and its redshift is just the 
non-relativistic Doppler shift z = v, leading to Hubble's law 

z(= v) = Horo (1) 

Hubble's law is well established because relative distances are easy to establish. All one has to do 
is find a 'standard candle', that is a type of object (say as star of a given type) of which all examples 
have practically the same luminosity. Then its apparent luminosity will vary with (distance) and 
so measure relative distances.^ On the other hand to fix Hq which is the constant of proportionality 
one has to know the luminosity of some object, which is much harder to do. Different estimates give 
Ho in the range 40 to lOOkmsec"^ Mpc~^, and it is usual to define a quantity h by 

iJo = lOO/ikmsec^^Mpc"^ (2) 

Thus, the redshift determination gives A < h < 1. As one might expect, Hq enters into many other 
equations besides Hubble's law, so the redshift determination is not the only method of determining 
Hq. I will suggest a 'best fit' value and error later on. 

Since H = a/a, the time taken for the universe to expand by an appreciable amount is of order 
the Hubble time 

H^^ = 9.7(3h-^ xlO^yr (3) 
During a Hubble time, light travels a distance of order the Hubble distance, 

Hf^^ ^ 2998h-^ Mpc (4) 

Hubble's law applies only to galaxies whose distance is much less than the Hubble distance, it is 
based on the non-relativistic Doppler shift which requires v 1. As we will discuss later, much more 
distant galaxies are also observed, the radius of the observable universe being of order the Hubble 
distance. Thus we are at the moment discussing only the nearby part of the observable universe. 
The universe has expanded by a negligible amount since light from the nearby universe was emitted, 
but by a significant amount since light from the distant universe was emitted. 

The Big Bang 

At the present epoch the Hubble time is of order 10^° yr. In what follows we shall extrapolate back 
to an era when the universe is very hot and dense, and the Hubble time is only a tiny fraction of 
a second. This era is popularly known as the Hot Big Bang, but one should not be misled by the 
language. Since the early universe is even more homogeneous and isotropic than the present one the 
expansion cannot be said to originate from a central point. It is certainly not an explosion, which 
by definition is driven by a pressure gradient. Thus the undoubted fact of the Big Bang does not 
explain the expansion of the universe; it has to be laid down at the very beginning. 

When is the beginning? Presumably it lies at the Planck epoch, when the Hubble time is of order 
the Planck time 

tpi = G^/^ ^ 5 39 X 10-44 (5) 

As we extrapolate back to this epoch, quantum gravity effects presumably invalidate the concept of 
time, which conveniently removes the need to discuss what is before the Big Bang! 

^To establish that an object is a standard candle one should ideally understand it theoretically, but in practice it 
is fairly safe to assume that a sufficiently distinctive object is standard. As a check one can compare two different 
candles, observing pairs of them which are known for to be in the same region of space because, for instance, they are 
in the same galaxy or star cluster. If they both have luminosity proportional to (distance)"'^, either they are both 
standard or they conspire to be non-standard in the same way. 



1.4 The matter content of the universe 



If gravity were negligible, a(t) would increase linearly with time and d would be constant. In fact, 
gravity slows down the expansion making d a decreasing function of time. To calculate the effect of 
gravity, consider a test particle of unit mass, on the surface of a comoving sphere (one expanding 
with the universe). If p is the mass density of the universe and r is the radius of the sphere, the 
mass inside it is (47r/3)r'^p and the potential energy of the particle is ~ {An / 3)r^ pG / r . The kinetic 
energy of the particle is and therefore 

P/2 - (47r/3)rVG/r = E (6) 

where E is the total energy of the particle. Writing r{t) = a{t)x and dividing through by a? /2 one 
arrives at the Friedmann equatiom 

H -—P-^ (7) 

where K = —2E. If A' < (energy E > 0) the expansion will continue indefinitely whereas if K > 
it will eventually give way to contraction leading to a 'Big Crunch'. The critical value separating 
these possibilities is K ^ 0, which corresponds to Friedmann equation 

= f^. (S) 
The corresponding mass density is called the critical density, and its present value is 

PcO = 1.88 X 10-29 /i2g(,j^-3 (g) 

= lO.S/i^GeVm-^ (10) 
= (3.0 X 10-3 eV)'^/i2 (11) 

It is convenient to define the density parameter by 

n =. p/p, (12) 

Since mass is conserved, p cx a^^. Putting this into the Friedman equation and remembering 
that H = d/a, we have an expression for dt/da which can be integrated to give t{a) and hence a(t). 
The result depends on K, or equivalently on the present density parameter Qq. As I shall discuss 
in Section 6, inflation strongly suggests the value Oq = 1, and almost everyone working in the field 
believes that this is the true value. ^ In that case Eq. (H) gives a =oc t^/^. Another case which is 
easy to solve is i^o ^ 1, which corresponds to negligible gravity and a cc t. 



Baryonic matter 

What is the observed value of flo? Consider first the contribution of ordinary matter ie nuclei 
and electrons. In the context of cosmology this is usually called baryonic matter since the baryons 
(nuclei) vastly outweigh the electrons. From the nucleosynthesis calculation we know that the baryon 
contribution to fio is [§) 

rJs (•013± .002)/i-2 < 09 (13) 

where I have used .4 < ft, to obtain the upper limit. Thus if fJo = 1 there exists non-baryonic dark 
matter, whose nature I discuss later. 

The luminous matter in the universe, consisting of stars and radiation-emitting gas, accounts for 
only fls ^ .01. Unless h is close to 1 (which we shall see in a moment is impossible with the favoured 
value total value CIq = 1), it follows from Eq. ( p^ that there is a lot of baryonic dark matter. If 
51o = 1 it constitutes on average a few percent of the total amount of dark matter. The fraction 
is presumably more or less the same throughout inter-galactic space. Within a given galaxy, one 

^It is often stated that Newtonian gravity cannot be used to discuss the expansion of the universe, and in particular 
that one cannot justify the neglect of matter outside the sphere in the above argument. This seems to me to be quite 
wrong. Drawing a much bigger sphere around the one that we considered, the matter within the bigger sphere can 
certainly be neglected. There remains the matter outside the sphere but it accelerates the small sphere and the test 
particle equally so it too can be ignored. This is exactly the same argument that justifies the use of Newtonian gravity 
for a galaxy or for the solar system, and it seems to be just as good in the cosmological case. 

^In this context Q includes all forms of energy density, but in the matter dominated era other forms of energy are 
supposed to be negligible except conceivably for that due to a cosmological constant as discussed later. 



might expect the baryons to be concentrated more in the central, luminous part than in the dark 
halo. The reason is that baryons (ordinary matter) can emit radiation whereas non-baryonic dark 
matter interacts too weakly to do so (or it would not be dark). In consequence baryons can lose 
more energy, allowing them to settle more deeply into the galaxy centre. The baryons in a galaxy 
might be in the form of a non-emitting gas, or they might be failed stars withQmass perhaps '--^ .001 
to .IMq, or dead stars (old white dwarfs, non-emitting neutron stars and the occasionally black 
hole arising from the collapse of an exceptionally massive star) with mass ^ Mq. These objects are 
called MACHOS (massive compact halo objects), because they will occur only in galaxy halos. (As 
far as we know, bound objects form only within galaxies so that the intergalactic baryons have to 
be in the form of non-emitting gas.) 

Detection of MACHOS in our galaxy has recently been claimed through their gravitational lensing 
of stars (microlensing) , a variation in brightness of the star occuring on a timescale of a few days 
as a MACHO crosses the line of sight If the observed MACHO density reflected the universal 
average baryon dark matter density, it would correspond to ^ -1, but since only the central part 
of the galaxy is probed this need not be the case, and even a much higher MACHO mass density 
would be quite consistent with the nucleosynthesis value — -05. 



Non-baryonic dark matter 

One can try to estimate the total amount of matter through its gravitational effect. The gravitational 
field in a bound system such as a galaxy or galaxy cluster can be deduced from the velocities of its 
components, as evidenced by the Doppler effect (the 'components' can be gas molecules, stars or 
whole galaxies) . One finds that each galaxy is surrounded by a dark halo accounting for most of its 
mass, the total galaxy contribution being CIq ~ .1. Galaxy clusters contain in addition inter-galactic 
gas, whose contribution to flo is not yet known. 

On larger scales, where the universe is almost homogeneous and isotropic, one can in a similar 
spirit observe the small departure from uniform expansion. This defines a 'peculiar velocity' field, 
which is usually called the bulk flow. If one knew the bulk flow v and the density perturbation Sp/p, 



one could deduce fio through the relation (||1|], cf. Eq. ( [llOD ) 



^ = -^(^0)-^ (14) 
3-Ho 3 p 

To estimate these quantities one has to rely on observations of the galaxies. Largely on the basis of 
numerical simulations of gravitational collapse, it is generally assumed that their motion accurately 
measures v, but that their density contrast is equal to a bias factor b times the underlying density 
contrast, with 6^1 within a factor of 2 or so (its value depends among other things on whether one is 
looking at optical or infrared galaxies). A recent study using this method ||] gives (ilo) *^ = 1.28t;5g67 
where bj is the bias factor for the infrared galaxies in the IRAS catalogue. Including non-linear effects 
strengthens this result slightly, so that for example if one assumes bj > .5 one deduces flo > -3 at 
95% confidence level. A smaller bias factor would certainly lead to problems with interpreting other 
data so it seems fair to say that the bulk flow indicates that fig > .1. In view of the nucleosynthesis 
limit ^Ib < .1 this means that non-baryonic dark matter seems to be needed. 

Further evidence about Qq comes from what one might call the 'standard model' of structure 
formation. This model, which has been the almost universal paradigm for many decades and is by 
far the most thoroughly explored and successful one, supposes that the structure originates as a 
Gaussian adiabatic density perturbation, whose spectrum at horizon entry is approximately scale 
independent. It requires that Qq is at the upper end of the range to 1, with — .1 definitely 
excluded. 



Summary 

From nucleosynthesis, baryonic matter contributes f^o — -01 to .09. The observed total Slo, especially 
if one accepts the standard model of structure formation, is higher, and the theoretically favoured 
value is r^o = 1. Thus non-baryonic dark matter is indicated. In addition, baryonic dark matter is 
indicated because luminous matter accounts for only fig .01. 

collapsing gas cloud has to have mass > Mq to achieve a high enough temperature for nuclear reactions to 
start, but on the other hand very light gas clouds are thought not to collapse at all. 



Value of fl at early times 



Even if is not equal to 1 at the present epoch, it quickly approaches 1 as we go back in time 
because the first term in Eq. (Q) is proportional to p oc a~^, and it therefore dominates the second 
term which is proportional to a^^. If fig 1 we can easily estimate the epoch before which ft is 
practically equal to 1 as follows. As long as O <C 1, gravity is negligible and therefore d is constant, 
leading to H (x a and fl cx p/H^ oc a~^. If follows that ~ 1 before the epoch 

a~f7o (15) 

Since flo > .1, we see that f2 is certainly close to 1 before a .1. 



The age of the universe 

One would like to define the age t of the universe at a given epoch as the time since a — (in other 
words, to set i = at a = 0). As we noted earlier the Planck scale puts an absolute limit on the 
earliest meaningful time, and more seriously the matter dominated era that we are treating here 
goes back only to a ~ 10"'^. So a practical definition of the age of the universe at a given epoch is 
simply the time since a was much less than the current value. 

The Friedmann equation allows one to calculate the present age of the universe in terms of flo. 
For Qq — 1 the behaviour a oc t'^^^ gives at any epoch 

and therefore 

to='^H-^ ^ 6.5x10^ h-^yi (17) 

The smallest conceivable value i^o — -1 gives 

to = .9i/(7^ = 8.8 X 10^h~^ yr (18) 

An upper limit on the age of the universe is provided by the age of the oldest stars (observed in 
globular clusters) which is bigger than 1.0 x 10^° years. With the favoured value fl^ — 1 this requires 
h < .65, whereas with fl^ = .1 the limit is h < .88. 



The measured value of the Hubble parameter 

We noted earlier that using Hubblc's law, different astronomers have produced estimates of h in the 
range A < h < 1.0. These Hubble law estimates are generally agreed to be very difficult, whereas 
the age limits that we have discussed are generally agreed to be sound. U Qq = 1, and if there 
is no cosmological constant (see below), one should therefore discard those Hubble law estimates 
which conflict with the age bound. Since both of these conditions are widely agreed to be likely, one 
concludes that the current best guess for the true value of the Hubble parameter is 

AO<h< .65 (19) 



1.5 High redshift observations 

So far we have discussed only the relatively nearby part of the universe, located at a distance small 
compared with the Hubble distance. The redshift Eq. (|l|) from galaxies in this region can be written 

z = da/a<l (20) 

where da is the change in the scale factor since the light was emitted. Thus, we see the nearby 
universe as it was in the very recent past, when the scale factor had almost its present value a = 1. 
Now we consider the distant universe, which is seen as it was in the distant past. 

In order to discuss the distant universe we need general relativity. The central idea of GR, and 
practically the only one that I shall use in these lectures, is that all measurements are to be made 
locally. Thus, one populates the universe with comoving observers. In the homogeneous, isotropic 
universe comoving observers are in freefall (their worldlines are geodesies) and as a result their clocks 
can be synchronized once and for all. In other words, they all agree on the age of the universe t. 
The hypersurfaces in spacetime with fixed t are orthogonal to the comoving worldlines, and will be 
called comoving hypersurfaces. 



The redshift measures the scale factor 



To calculate the redshift from a distant object, emitting its light when a <C 1, we can add up the 
redshifts seen by a sequence of comoving observers (observers moving with the expansion of the 
universe), each separated by a distance <C H^^. In the small region around each one Eq. ( po| ) still 
applies except that it is evaluated at the time t when the radiation passes by. Using H = d/a and 
integrating gives 

! + . = (21) 

A 

The redshift measures directly the scale factor of the universe at the time of emission of the radiation. 
The size of the observable universe 

We would like to know the present distance of an object observed with redshift z, defined as the 
sum of the present distances between the comoving observers along its line of sight. The distance 
dr from one comoving observer to the next is proportional to a(t) and so can be written dr — adx 
where dx is independent of time. Thus each observer has a time independent comoving coordinate 
X, related to his distance from us by r = ax. Let us work out the trajectory r{t) = x(t)a(t) of a 
photon moving towards us. During time dt it passes between observers with coordinate distance 
x + dx and who are separated by distance adx — dt {ie thes e observers measure the photon speed 
to be 1). The comoving distance of the source from us is therefore x = dt/a, where the subscript 
denotes the present and the subscript e the time of emission. Its present distance is therefore 

r(to,te)-ao^ dt/a (22) 

(It is convenient to display the present value ao = 1 of the scale factor in this context.) Given the 
time dependence of a, this relation gives the present distance of an object observed with redshift z. 

In the limit z ^ 1 it gives the present distance that we can see in principle, corresponding to 
the limit where the object emits its light at t = 0. This distance is called the particle horizon and 
is given by 

/■*" dt 

rp.h.{to) = ao / — (23) 
Jo ^ 

Assuming that JIq = 1 one has a oc t'^^^ and therefore ?'p.h.(io) = '^Hq^. 

There is nothing special about the present, so we can define the particle horizon i?p.h.(^) ^iny 
epoch, 



, , r dt , , r I da 
Jo a Jq aH a 



(24) 



It gives the size of the biggest causally connected region, across which a signal with the speed of 
light can travel since t = 0. 

As we noted earlier, the mathematical epoch a = t = has to be replaced in practice by an 
epoch when a and t are much less than their current values. Thus, rp±,(t) can really only be defined 
as r{t,te) where te is some early epoch. As long as gravity is attractive, aH = a increases as we 
go back in time, and from Eq. ( [2^ ) r{t,te) converges as te 0, and is typically of order H'^{t) as 
we just found for our particular case {Qq — 1 and matter dominating). But if the universe begins 
with an inflationary era, where gravity is by definition repulsive, r(i, t^) diverges. If we push t^, back 
to the beginning of infiation, it is much bigger than the Hubble distance H~^(t), so a region much 
bigger than the Hubble distance is causally connected.^ But to achieve this causal connection at a 
given epoch, one has to work with signals emitted at a much earlier epoch, and what usually matters 
physically is the size of a region that can be causally connected by signals emitted in the last Hubble 
time or so, namely the Hubble distance. In the context of the early universe, the Hubble distance 
is usually referred to as 'the horizon', and has largely displaced the 'particle horizon' as the sort of 
horizon that cosmologists are interested in. 

^At all epochs t up to the present, on the usual assumption that the observable universe is inside the Hubble 
distance at the beginning of inflation. 



Non-Euclidean geometry? 

To go further in interpreting high redshift observations, one needs to know the spatial geometry, ie 
the geometry of the fixed t hypersurfaces. Einstein's field equation gives the spatial line element as 

dr^ 

= E^TTT + (d^' + sin' Od^^ ) (25) 

where K is the quantity that appears in the Friedmann equation. The distance element dl is that 
measured by a local comoving observer. 

A surface with fixed coordinate r clearly has the geometry of a sphere in Euclidean space with 
radius r, and in particular its area is Anr^. However the radial distance between adjacent spheres 
is (1 — Kr'^ /a^)~^/'^dr, which is equal to dr only if K = 0, corresponding to f2 = 1. Thus space 
is non-Euclidean ('curved') ifflj^l. If < 1 space is infinite just as in the Euclidean case, but 
if f2 > 1 it has the remarkable property of being finite but unbounded. In particular, a surface 
made out of geodesies, such as the 'equatorial plane' 4> = 7r/2, has the geometry of a sphere whose 
curvature is — Kja^. The quantity R is the radius that the sphere would have if it lived in 
Euclidean space, or equivalently it is defined by saying that the area of the sphere is AirR^. The 
total volume of space, obtained by integrating the volume Attt'^ {1 — Kr^ / a'^)~^^'^ dr between adjacent 
spheres, is Idir'^R'^. 

One calls K/a^ the 'curvature of space' whether it is positive (f2 > 1) or negative (fi < 1). The 
departure from Euclidean geometry becomes significant when > \a'^/K\, and one refers to the 
distance \a^/K\^^^ as the curvature scale. If ~ 1 the curvature scale is much bigger than the 
Hubble distance H^^, if |1 — f7| <C 1 it is of order H^^, and if 3> 1 it is much less than H^^. 



High redshift observations 

Having displayed the basic tools I will briefly say what we learn from observation of the distant 
universe. The most important finding is that the universe is definitely evolving. The most dramatic 
case is that of quasars (active galactic nuclei), whose abundance per comoving volume peaks at z ~ 3 
or so. Neither quasars nor any other objects are observed at z > 5. Ordinary galaxies as well as 
clusters are observed out to a redshift of order 1 to 2, and they too show signs of evolution. (The 
mere fact that galaxies have size of order one tenth of their spacing and do not expand with the 
universe means, of course, that they cannot have existed before z ~ 10.) 

Ideally, high redshift observations plus an understanding of galactic evolution would give us 
information the value of f2o- In practice the information obtained in this way does not add anything 
to our knowledge as yet, except for the issue of a cosmological constant that we look at next. 



The Robertson- Walker metric 

Before leaving the description of the large scale universe I should write down the famous Robertson- 
Walker spacetime metric, 



dx^ 



+ x^{de^ +sm^ 0d(j3^) 



(26) 



It expresses the fact that the universe is homogeneous and isotropic, with spatial geometry given by 
Eq. (p5|). The spacetime metric is the traditional starting point for cosmology but in these lectures 
I will consistently avoid it, in favour of the physical notion of comoving worldlines. 



1.6 A cosmological constant? 

The discussion so far assumes that a) gravity is described by Einstein's field equation, and b) the only 
significant source of gravity is the matter. People have considered the possibility of alter assumption 
a) by adding a cosmological constant A to the lagrangian of Einstein gravity. This is equivalent to to 
keeping Einstein gravity, but modifying assumption b) by invoking an additional source of gravity 
in the form of a fluid with pressure pA and energy density pA related by p = —p. This second 
viewpoint is the most useful and I adopt it here. From Eq. (^) below it then follows that pA is time 
independent, so we are considering the possibility that 



p{t) = p„i{t) + PA 



(27) 



The corresponding present density parameter is 



Qo = ^rn + (28) 

If we are ever faced with a measured value flm < 1, we will probably want to invoke a cosmological 
constant to keep the total value fio = 1 which inflation suggests. A cosmological constant is, however 
unattractive from a particle physics viewpoint because the corresponding tiny vacuum energy density 
{.OOleV)'^ (Eq. (pl|)) is difficult to understand. What about observation? Assuming that the total 
Oo = 1, the present the situation is as follows 0, 1. 

On scales much less than the Hubble distance the effect of a cosmological constant is to add a 
repulsive term to Newtonian gravity, but the extra term is insignificant. 

Going on to observations of the distant universe, gravitational lensing statistics give a remarkably 
good limit fl^ < .6 1^, and structure formation gives a similar but less reliable limit. Thus > -4, 
which means that one cannot use a cosmological constant to avoid the need for non-baryonic dark 
matter. 

Coming to the age of the universe, Eq. ( ^ ) shows that ft^ corresponds to repulsive gravity, 
a > 0, so the age of the universe will be an increasing function of £7 a- In fact, the upper limit on 
h increased from .65 to .88 as ftj^ increases from to .6. Thus a high value of h measured by, say, 
Hubble's law would be evidence in favour of a cosmological constant, and conversely a confirmation 
of the expected low value would be evidence against it. 

In summary, a cosmological constant is theoretically objectionable, and is not so far required by 
observation. 



2 The Radiation Dominated Era 
2.1 Overview 

An epoch in the early universe is conveniently specified by giving the temperature T in eV, MeV 
or GeV (as usual in cosmology I set fcs = !)■ The epoch from 1 MeV to leV, which can be taken 
to mark the end of the early universe, is well understood, and is highly constrained because it 
includes the epoch of nucleosynthesis. At earlier epochs one encounters a lack of knowledge about 
the nature of the quark-hadron and electroweak phase transitions, and eventually a basic ignorance 
of the fundamental interactions and the initial conditions. Later one enters the matter dominated 
era, when structure formation occurs and again our knowledge becomes less secure. The various 
epochs are summarised in Table 1. 



2.2 The cosmic microwave background (cmb) 

The most direct evidence for a hot early universe comes from the cosmic microwave background 
(cmb). It is extremely isotropic, with a very accurate blackbody distribution corresponding to 
temperature T = (2.736 ± .017)° K. 

To understand what the cmb implies, recall that the blackbody distribution specifics the photon 
occupation number of each momentum state as 



/(p) = 2[e(^'/^)-l]-i 



(29) 



The factor 2 counts the number of spin states, and p is the momentum. Since there are (27r)-3^V^ 
states in a given volume of phase space, the number density n and the energy density p of the photons 
are given by 



P = 



(27r)3 
2 

(2^ 



Ef{p)4:Trp'^dp 



(30) 
(31) 



with E = p. Evaluating the integrals gives 



Pi 



2(7rV30)r^ 
2(C(3)/7r2)r3 



(32) 
(33) 



TEMPERATURE 


EVENT 


CONSEQUENCES 


mpi = lO^'-'GeV 


Inflation begins? 


Immediate collapse avoided 


lO^^-lO^^ GeV? 


Cosmological scales 


Field fluctuations 




leave the horizon 


become classical 


lO^^-lO^^ GeV? 


Inflation ends 


Particles created 


10^^ GeV-1 TeV?? 


Particles thermalize 


Radiation domination begins 


lTeV-100 GeV 


Electroweak phase 


Particles acquire mass 




transition 


Baryogenesis occurs? 


100 MeV 


Chiral and quark/hadron 


Pions acquire mass and 




phase transitions 


quarks bind into hadrons 


IMeV 


Neutrinos decouple 


Neutrons start to decay 


IMeV 


ec annihilation 


Afterwards (T^/T^) = (ll/4)i/3 


.IMeV 


Nucleosynthesis 


Most neutrons bind into "^He 


1 eV 


Matter domination 


ihe cold dark matter density 




begins 


contrast starts to grow 


.3eV 


Atoms form 


Photons decouple and the baryonic 
matter density contrast grows 


~ 10"^ eV? 


First galaxies form 




2.4 X 10"* eV 


The present epoch 





Table 1: The early universe. The left hand column defines the epoch. It gives the temperature, 
which is roughly equal to p^^^ where p is the energy density (in the first three rows the universe is 
not actually thermalised so the number in the left hand column is defined as p^''*). During radiation 
domination the time since the big bang is roughly 1 MeV/T seconds, and during both epochs the 
scale factor is roughly a ~ Tq/T where Tq ~ 10"'' eV is the present temperature of the cosmic 
microwave background. 

where C(3) — 1.2021. The photon energy density is negligible compared with the mass density of 
the matter, corresponding to present density parameter 

= p^/pco = 2.51 X lO-'^/i"^ (34) 

Dividing n by the baryon number density ns — 11.2h? m^^ils, 

ri= — = 2.68 X IQ-^h^VlB = (3.4 ± .6) x 10"^" (35) 

(The last equality is the nucleosynthesis result Eq. (|l^).) Thus there are many photons per baryon. 

Now let us evolve the blackbody distribution back in time, assuming that the number of photons 
in a comoving volume is conserved (as we shall see this is correct back to the epoch of electron- 
positron annihilation). Because of the redshift, the momentum of each photon goes like a~^, which 
is equivalent to saying that the blackbody distribution is preserved, with T cx a~^. Since T cx a~^, 
p^ oc a^^, so Q,^i/Q,m oc a~^. We conclude that before the epoch a^^ ^ 10"* the energy density of the 
cmb photons dominates that of the matter. 

2.3 Physics in the radiation dominated era 

Motivated by the observed cmb, one makes the hypothesis that the early universe is a hot ideal gas. 
The hypothesis is found to be consistent in the context of gauge theories like the Standard Model. 
It leads to the prediction that the nucleosynthesis of ^He and a few other light nuclei occurs in the 
early universe, and the calculated primordial abundances are in striking agreement with observation. 

Adiabatic expansion 

The pressure of a gas is given hy p — u^p/3, where is the mean square particle velocity and p is 
the energy density. It follows that 



p — p/3 relativistic (36) 

j3 ~ nonrelativistic (37) 



As the universe expands, the change in the energy E = a'^p in a comoving volume a"^ is given by 
dE = ~pd{a^) so that 

a^--3(p + p) (38) 

For a nonrelativistic gas this gives p oc a~^, but for a relativistic gas it gives 

p cx a-^ (39) 

The extra factor arises from the redshift of the particle energies between collisions, each collision 
conserving energy. 

The homogeneity and isotropy of the early universe implies that its expansion is adiabatic (no 
heat flow), so that entropy is conserved. Keeping only relativistic particles in the gbb, the entropy 
density is 

P + P 4p 4 7r2 3 
^ = ^-3T = 3 30^*^ 

Conservation of the entropy a^s in a comoving volume therefore gives the relation between the 
temperature T and the scale factor a 

aToig^^^^ (41) 

For rough estimates one can ignore the time dependence of to deduce that T cx a^^, and can 
also set s ~ n to deduce that particle number is roughly conserved. In particular, Eq. ( ^Sf ) holds at 
least roughly in the early universe, 

ne/n^ - 10"^° (42) 

For a detailed calculation one has to take into account that aT increases abruptly whenever the 
temperature falls below the mass of a particle species. Physically, the relative increase in temperature 
comes from radiation emitted when the species annihilates and/or decays. 

The Friedmann equation 

According to general relativity, the Friedmann equation Eq. (^ that we derived from Newtonian 
gravity remains valid, if p denotes the energy density instead of just the mass density. One can 
differentiate it using Eq. (pq), which is equivalent to 



H = -H^-—{p + 3p) (44) 



p=-3H{p + p) (43) 

The result is 

H = - H^ - - 

3 
or 

d AttG , , , , 

- = ^P + 3p 45 

a 3 

This last expression reduces to the Newtonian result in the limit p <C p, if p is identified with the 
mass density. But in general both energy and pressure are sources of gravity. 

Returning to the Friedmann equation, the density parameter f2 is still defined by Eq. (p^). 
Setting it equal to 1, corresponding to = 0, and remembering that for radiation p oc a~^, one 
finds 

a oc (46) 

corresponding to 

(47) 

Thermal equiUbrium 

In the early universe collision and decay processes are continually creating and destroying particles. 
Let us suppose that it is in thermal equilibrium, which means that each process is taking place at 
the same rate as its inverse. Then the the number of particles of a given species, per momentum 
state, is given by 

/(p)=.9[e(^-^)/^±l]-i (48) 



In this expression, g is the number of spin states, p is the momentum, and E is the energy, given 
in terms of the mass m hy E ^ + to2_ The sign is + for fermions (ensuring that there is 



at most one particle per quantum state), and — for bosons. The quantity fi (which depends on 
temperature) is called the chemical potential of the species. Given the chemical potential, this 
expression determines the number density n and energy density p of the species through Eqs. ( |30| ) 
and @. 

The chemical potential is conserved in every collision, so that if a process AB — > CD occurs 
then Ha + fJ-B — fJ-c + fJ-D- In the early universe all known particle species are freely created and 
destroyed (provided that their mass does not exceed the typical collision energy, which we shall see 
in a moment is of order T). The only significant restriction is that each collision must respect the 
conservation of any 'charges' that the particles carry. In the context of the Standard Model there 
are five such charges, namely electric charge, baryon number and the three lepton numbers]^ Since 
the photon carries none of these 'charges', a single photon can be created through processes of the 
form A — » A7, which implies that fi = for the photon, leading to the blackbody distribution. The 
same goes for any particle which is its own antiparticle, so it too has /i = 0. If the antiparticle is 
distinct one can create a particle-antiparticle pair through some process of the form A —^ A+pair, 
which implies that the particle and its antiparticle have opposite fi. As a result, /i vanishes if the 
number densities n and n of particle and antiparticle are equal, and otherwise is determined by the 
imbalance n — n. 

More generally, one can show that until weak interaction equilibrium fails at T 1 MeV, enough 
reactions occur to determine all of the chemical potentials in terms of the densities of the five 
conserved 'charges'. Furthermore, all of the chemical potentials are zero if the 'charges' are all zero. 
In that case one has what I shall call the generalized blackbody distribution or gbb 

f{p)=9[e''^^±l]-' (49) 

Putting this into Eqs. ( ^0| ) and (|l|) gives n and p. If the temperature is well above the mass of the 
species in question, m <^T, then one is in relativistic regime and it is a good approximation to set 
E ^p. Using Eqs. ^) and (|l]) gives 

_ r (7rV30)gr4 bosons 

^ ~ \ (7/8)(^V30)gr4 fermions ^^^> 

(C(3)/7r2)gr3 bosons 
(3/4)(C(3)/7r2)gr3 fermions ^ ' 

Thus according to the gbb, each relativistic species contributes ^ to p and ~ to n. (Note that 
the typical energy E ~ p/n is of order T, so this is indeed the relativistic regime.) As T falls below 
m one moves into the nonrelativistic regime, and according to the gbb p and n fall rapidly (like 
g-m/T^_ The physical reason is that the typical energy E available in a collision is now insufficient 
to create the species. 

Is the gbb valid in the early universe? The charge density of the universe is certainly zero to very 
high accuracy, or the expansion of the universe would be governed by electrical repulsion instead of 
gravity. The baryon number is not zero, but it is small in the sense that 

T] = nB/n^ (52) 

The three lepton number densities cannot be measured directly, but let us suppose that they too 
are small in the same sense. It can be shown that if this is so, the gbb is valid to high accuracy for 
all relativistic species, the reason being that for each species one then has \n — n\ T^. From 

now on I assume that the gbb is valid for all relativistic species in equilibrium, and the success of 
the nucleosynthesis calculation will support that assumption. 

As the temperature falls below the mass to of a given species, particle-antiparticle pairs at first 
rapidly annihilate as dictated by the gbb. Then the small particle-antiparticle imbalance becomes 
significant, and annihilation soon stops because only particles (or only antiparticles) survive. Even 
if they do not decay, the surviving particles give a negligible contribution to p during radiation 
domination, and an even more negligible contribution to n. 

As long as a non-relativistic particle species is in equilibrium, it follows from Eqs. ( ^8|) and ( |30| ) 
that its number density is 

" = (53) 



®The non-perturbative violation of the last four is not relevant in this context. 



We shall use this equation to determine the relative numbers of different species of nuclei in thermal 
equilibrium, and the degree of ionization of a gas in thermal equilibrium. 

The upshot of this discussion is that unless one happens to be at an epoch when T is close to 
one of the particle masses, p and n come just from the relativistic species. In particular, 

p= (7rV30)5*(T)r4 (54) 

where the effective number of degrees of freedom g, is given by 

9*iT)^ 9 + 1 E 5 (55) 

bosons fcrmions 

where the sums go over particle species with m < T. According to the Standard Model it is of 
order 10 for 1 MeV < T < 10 MeV and of order 100 for T > 1 GeV. Extensions of the Standard 
Model, notably supersymmetry, introduce new particles which can increase g» by a factor of a few 
at T > 100 GeV, but usually leave it unchanged at lower energies. 

Eq. (^) applies only if all relativistic species are in thermal equilibrium. This is true in the 
Standard Model (for T > 1 MeV), but not necessarily in extensions of it. 

The timescale 

Using Eqs. (^J), (||) and (p7|), the timescale during radiation domination is 

2.4 The Standard Model of the early universe 
Interaction rates and the ideal gas condition 

Our basic assumptions are that the early universe is an ideal gas, and that it is in thermal equilibrium. 
The first requires that the interactions are not too strong, but the second requires that they are 
strong enough to make the relevant reaction rates (per particle) be bigger than the rate of expansion 
H of the universe. Which reactions are 'relevant' depends on what one supposes about the chemical 
potentials; each reaction that is in equilibrium implies a relation between the chemical potentials (if 
it creates or destroys particle species) and there have to be enough relations to ensure the postulated 
conditions on the chemical potential. In particular, if we want the gbb to hold for relativistic particle 
species, enough interactions have to be in thermal equilibrium to ensure that all of the chemical 
potentials vanish (on the assumption that the 'charge' densities all vanish). In the context of the 
standard model, one can show that this is indeed ensured if all processes of the form AB CD 
going by the exchange of a single gauge boson are in equilibrium. Let us check that such reactions 
are indeed in thermal equilibrium, and that the interactions are yet weak enough to ensure the ideal 
gas condition. 

Suppose first that the gauge boson is massless; this covers the case of the electromagnetic interac- 
tion before photon atoms form (photon decoupling) the strong interaction before hadrons form (the 
QCD phase transition) and the weak interaction before the W and Z acquire mass (the electroweak 
phase transition) . For an isolated collision the cross section a is infinite because the interaction has 
infinite range, but in a gas there is a cutoff at the particle spacing (Debye shielding). Using 

n ^ one finds that a ~ /T'^, where a is the electromagnetic, weak or strong coupling strength 
as appropriate. If at least one of the initial particles is relativistic, the rate per particle is F ~ na 
(the formula is exact if one particle is at rest, and a factor 2 too small for a head-on relativistic 
collision.) Using Eqs. (^) and (|5^) this gives T/H ~ .la^nipi/T. Taking the high energy estimate 
a '--^ .1, this is indeed ^ 1 unless T > 10^^ GeV. Thus the interaction is strong enough to maintain 
thermal equilibrium, except at these very high temperatures. What about the ideal gas condition? 
We can take it to be the requirement that the mean free path r, determined for relativistic particles 
by ran ^ 1, is big compared with the particle spacing n~^^^. Thus we require n^^^a <C 1 which 
is equivalent to a ^ 1. This is satisfied, except for the strong interaction just before the QCD 
transition. 

Now consider the weak interaction after the electroweak phase transition. The cross section is 
a ^ a^T"^ /m^r, which gives T/H ^ mpia^T^ /m^^. Taking the weak interaction coupling as a ^ .01 
this is equal to 1 when 

T - a-2/3(„^4^/TOp;)l/3 _ I to 10 MeV (57) 



Weak interactions are in equilibrium only above this temperature. Since the cross section is a factor 
(T/mw)^ smaller than in the previous case, the ideal gas condition is amply satisfied for the weak 
interaction. 

2.5 The history of the early universe 
The electroweak phase transition 

The electroweak symmetry of the standard model is restored above some critical temperature Tc ~ 
100 GeV. The Higgs field is then zero, corresponding to a false vacuum with energy density ~ 
(100 GeV)*, and all particle masses vanish. As the temperature falls through Tc the Higgs field 
rapidly settles to its true vacuum value and particles acquire mass. 

The consequences of a phase transition can be very different, depending on whether it is of 
first order or second order. In the first case the transition is relatively violent because bubbles of 
the new phase are produced. In the second case it is more gently, the new phase being produced 
smoothly everywhere (except near topological defects, which are not produced at the electroweak 
phase transition). 

It is not clear if the electroweak transition is of first or second order, and the answer in any 
case should probably take into account an extension of the Standard Model. If it is of first order it 
will violate baryon- and lepton number conservation because of non-perturbative effects occurring 
in bubble walls. As Shaposhnikov explains in his lectures, such a mechanism might be responsible 
for generating the observed ratio rib/n-y 10~^°, either within the standard model or within some 
extension of it. 

The quark-hadron transition 

As the temperature falls through Aqcd ~ 100 MeV, chiral symmetry is spontaneously broken which, 
among other things, gives mass to the pion. More importantly, the quark-hadron phase transition 
occurs, which binds quarks and gluons into hadrons as their typical spacing ^ T^^ becomes bigger 
than AqJjd ~ 1 fm. . . ^ . 

The quark-hadron phase transition H is thought to be a second order one. If it turned out to be 
first order though, the bubbles could lead to dramatic effects, perhaps interfering with nucleosyn- 
thesis or even creating black holes accounting for the dark matter. The mass of such black holes 
might be of order IMq corresponding to Eq. (|7|) below, or much smaller. 

Neutrino decoupling and ee annihilation 

At 1 MeV <T< 100 MeV, the only particles present (according to the Standard Model) are photons, 
free protons and neutrons, electrons, positrons, and the three neutrino species plus their antiparticles. 
Other elementary particles as well as hadrons are absent because they are unstable, and are too heavy 
to be created in a typical collision (to 3> T). Nuclei, not to mention atoms and molecules, are absent 
because their binding energy is much less than T, so that collisions would destroy them. 

The protons and neutrons turn into each other through reactions like ep t'eU, en ^ D^P and 
n pei^e- As a result the ratio n/p of their number densities is determined by thermal equilibrium. 
The chemical potentials satisfy pLp — ^in = y^v^ — Me, with /x^^ and /ie both <C T (on the assumption 
that the lepton number densities are sufficiently small that the e and Vf. satisfy the gbb to good 
accuracy). From Eq. ( |5^ ) with |/Zp — m„| <C T, it follows that the neutron to proton ratio is 

n/p = e-^"/^ (58) 

where Am = to„ — rup = 1.3 MeV. 

When the temperature falls through 1 MeV, three essentially unrelated things happen to occur 
at about the same time. 

• Neutrinos decouple because of Eq. ( p7| ) 

• Electron-positron annihilation occurs because TOg ^ 1 MeV 

• Neutron creation stops because of Eq. ( |57| ) .[| 

^It would stop anyway, because to create a neutron requires a collision energy bigger than m„ — rrip = 1.3 MeV. 



A detailed calculation shows that the events actually occur in the indicated order; let us discuss 
them in turn. 

After decoupling the neutrinos travel freely, maintaining their gbb with temperature Tj, cx a~^. 
The photon temperature, though, is relatively raised by ee annihilation. Before annihilation there 
are electrons, positrons and photons in equilibrium giving g^, = 11/2, whereas after there are only 
photons giving = 2. From Eq. ( p] ) it follows that after annihilation, = (ll/4)^/^T^. Re- 
membering that there are two spin states for the photon and one for each massless neutrino or 



antineutrino species, Eq. (51) gives 

= (3/ll)n^ (59) 



for each neutrino species. If all three neutrino species are massless, Eqs. ( pS] ) and (51) give the 
neutrino contribution to Qq, 

= (21/8)(4/ll)''/^r2.^ = .681f7.^ (60) 
Using Eq. (^) the epoch of matter-radiation equality is therefore 

Zaq = 2.37 X lO^/l^ (61) 

After neutrons cease to be created, the np ratio is frozen in except for the slow decrease caused 
by neutron decay. Its initial value is n/p ~ 1/6, and by the time that nucleosynthesis occurs at 
T ~ .1 MeV it has fallen to 

n/p -1/7 (62) 

Nucleosynthesis 

A nucleus is defined by the number Z of protons, and the number A — Z oi neutrons. If strong 
nuclear reactions occur at a sufficiently high rate, while the weak interaction changing protons to 
neutrons is negligible, the chemical potential of a given nucleus satisfies fi — Z^ip + {A — Z)fin- 
By virtue of this relation, the /i's cancel in the ratio n-^n^~'^ /n given by Eq. ([5^), leading to the 
following number density for a given nucleus 

/ o_ \3(A-l)/2 

n = 5^3/^2-^ (^-^j n>^^e^/^ (63) 

where B is its binding energy. Together with the np ratio and the total baryon density ub this 



determines all of the nuclear densities. Using ub = rjUj together with Eq. (33) for the photon 
number density n-y, one can calculate the mass fraction in a given nucleus. 



X = 5[C(3)^-i7r(i-4)/22(3A-5)/2]^5/2^ 

X X^^Xt'i^] e^/^ (64) 



If the baryon to photon ratio rj were of order 1, this expression would say that the nuclei form {X 
becomes of order 1) at roughly T ^ B, which is the epoch after which a typical photon has too little 
energy to destroy the nucleus. But because there are ri~^ ^ 10^*^ photons per baryon, the nuclei 
do not form until T is well below B. The reason is that until that happens there are still a lot of 
exceptionally energetic photons around with energy of order B. 

To determine whether a given set of nuclei will actually be in thermal equilibrium at a given epoch 
in the early universe, one has to take the rates of the relevant nuclear reactions from accelerator 
experiments. (They cannot all be calculated accurately from first principles, and are certainly not 
amenable to the order of magnitude arguments that work for elementary particles.) One finds that at 
T > .3 MeV the lightest few nuclei are indeed in thermal equilibrium, but that their number densities 
are negligible. In the range .IMeV <T< .3 MeV, thermal equilibrium predicts that practically all 
of the neutrons are bound into ^He, but when T first enters this range thermal equilibrium actually 
fails badly because the reactions forming "'He do not occur fast enough to keep up with the 'demand' 
for "^He (the deuterium bottleneck). As a result most of the neutrons bind into ^He at T ~ .1 MeV, 
rather than at T ~ .3 MeV as nuclear equilibrium would predict. Assuming that all of the neutrons 
bind into "^He, its abundance by weight is 



(65) 



Figure 1: Nucleosynthesis predictions and observations, taken from Q. The upper section gives 
the ''He abundance by weight, and the lower section gives the other abundances by number. The 
horizontal lines indicate the observed values; from top down the upper and lower limits for ^He, the 
upper limit for D+'^He, the lower limit for ■^He and the upper and lower limits for ^Li. The horizontal 
axis is lO'^Ty, and consistency with all observations requires that it be in the shaded region. 



Taking n/p ~ 1/7 gives X ~ 22%. 

The assumption of thermal equilibrium gives a reasonable estimate of the ^Hc abundance, but an 
accurate calculation of it and even an estimate of the other abundances requires an out of equilibrium 
calculation. This has to be done on a computer, taking the rates for the various nuclear reactions 
from laboratory experiments. The latest results from such calculations 0] are shown in Figure 1 as 
a function of 77. 

The primordial abundances of ''He and other light nuclei can be measured by a variety of tech- 
niques, and the predicted abundances agree provided that 77 is in the interval indicated, which 
corresponds to the result already quoted in Eq. (|3^). The agreement is extremely impressive, and 
constitutes one of the strongest reasons for believing the hypothesis that went into the calculation. It 
also limits greatly ones freedom to alter them; in particular it severely constrains possible extensions 
of the Standard Model. 



Photon decoupling 

After nucleosynthesis there are photons, protons, ''He nuclei and electrons in thermal equilibrium. 
When the temperature falls sufficiently, most of the electrons bind into hydrogen atoms, and one 
can use the thermal equilibrium fraction of hydrogen to determine when this happens. The atoms 
form through the process p + e — > H + 7, so that /ip + /Xg = ^lh- The chemical potentials cancel in 



the ratio npUe/nH given by Eq. (53), leading to 

-3/2 

27r", 

Here B = mp + rrie — ma — 13.6 eV is the binding energy of hydrogen, and I have used gp ~ ge ~ 2 
and gn — 4. Ignoring the ''He for simplicity, rig — Up and ub = np + nn, which with Eq. 



uh = ripUe — — e ' (66) 



determines the fractional ionization X — Up/uB, 



4V2C(3) /^y/^^,/, (^^^ 



This is cahed the Saha equation. As with nuclei, very few atoms form until T is well below 
because there are so many photons per baryon. In fact, one finds that T ~ .3eV when X has fallen 
to 10%. The corresponding redshift is z ~ 1300. Soon after this, at z ~ 1100 to 1200, the photons 
decouple. 



2.6 Beyond the Standard Model 

There are several theoretical reasons for wanting to extend the Standard Model. 

• One-loop contributions to the Higgs mass have to be cut off at the scale 100 GeV to 1 TeV, 
or they will be too big. To impose the cutoff one can either make the Higgs composite, as 
in technicolour type models, or invoke supersymmetry. Composite models are becoming in- 
creasingly unviable in the light of the ever more impressive agreement between the Standard 
Model and collider experiments, notably at CERN, and I will not consider them here. Super- 
symmetry requires that each known particle species has an as yet undiscovered partner. The 
lightest supersymmetric partner is typically weakly interacting and stable, and is therefore a 
dark matter candidate. 

• Unless one of the quark masses vanishes, a global U{1) symmetry called Peccei-Quinn sym- 
metry is presumably needed to explain why the strong interaction is CP invariant. To avoid 
conflict with accelerator physics and astrophysics, the symmetry is broken at a very high energy 
scale > 10^° GeV. 

• Going to a still higher scale, one would like to make contact with gravity, presumably involving 
superstring degrees of freedom around the Planck scale ^ 10^^ GeV. 

• One may wish to invoke a GUT. Superstring theory can apparently relate couplings etc without 
one, but on the other hand the ratio of the strong, weak and electromagnetic interaction 
strengths predicted by supersymmetric GUTS agrees with experiment amazingly well, with a 
unification scale ~ 10^^ GeV. 

An observational reason for wanting to extend the Standard Model might be the desire to introduce 
neutrino masses to solve the solar neutrino problem as well as a problem with atmospheric neutrinos. 

Because the cosmological predictions of the Standard Model agree so well with observation, many 
otherwise reasonable extensions of the Standard Model can be ruled out. Particularly sensitive are 
nucleosynthesis, and the absence of significant corrections to the shape of the cmb spectrum. 

From the many things that could be discussed, I choose three of the most topical. 



Neutrino masses 



A massive neutrino species could be a dark matter candidate. From Eqs. (11), (|33| ) and (^9|), its 
contribution to the density parameter is 



A mass of order 10 eV is therefore required. Experiments observing /3-decay practically rule out such 
a mass for the Vg, but it is allowed for the or Vr- 

A firmer reason for wanting some neutrino species to have mass comes from the solar neutrino 
problem, whose favoured solution requires either the or the i^r to have a mass of about 3 x 10^^ eV 
(with the Ve mass much smaller). 

A natural mechanism for generating neutrino masses is the see-saw mechanism, which invokes a 
right-handed neutrino with very large mass M , related to the lepton mass m and neutrino mass 
of each generation by 

= m^/M (69) 

Taking, without very good justification, M to be the same for all three generations, the choice 
M ~ 10^^ GeV gives m^^ ~ 10 eV and m^^ ~ 3 x 10""^ eV, so that the first neutrino can be the dark 
matter and the second can solve the solar neutrino problem. 



Neutrino dark matter is called hot, because it remains relativistic until the epoch when most 
cosmological scales of interest have entered the horizon (come into causal contact), and therefore 
cannot initially undergo gravitational collapse. 



Cold and warm dark matter 

Cold dark matter is by definition non-relativistic when all cosmological interesting scales enter the 
horizon. The lightest supersymmetric particle is a cold dark matter candidate. Its interaction 
strength is similar to that of the neutrino, but because it is much heavier it will be non-relativistic 
when it falls out of equilibrium. Its contribution to D,q turns out to be 

"--(liv)'""^ (™> 

Since one needs m ~ 100 GeV to 1 TeV for supersymmetry to work, this is automatically of order l! 

Warm dark matter by definition remains relativistic until a cosmologically interesting epoch, 
which is however significantly earlier than the epoch for neutrino dark matter. Candidates include 
a right handed neutrino, a majoran, or a particle with interaction strength very much less than that 
of neutrinos with a mass of order 1 keV. 

The dark matter candidates discussed so far are known collectively as WIMPS (weakly interacting 
massive particles). One can hope to detect a WIMP in the laboratory by searching for its rare 
interactions, the detectors being of the same type as the ones that see solar neutrinos or do neutrino 
astronomy (so far only with supernova 1987A), and place limits on proton decay. One might also 
hope to detect WIMPS in the galaxy by looking for photons produced by their annihilation. Finally 
one might be able to create them at accelerators. 

A quite different dark matter candidate is the axion, which is never in thermal equilibrium and 
would be cold dark matter. The axion is the Goldstone boson of Peccci-Quinn symmetry, which is 
the best bet for ensuring the CP invariance of the strong interaction (except perhaps for making 
one of the quark masses zero) . Various possibilities exist for axion cosmology in different regimes of 
parameter space |p^ , but typically axions are emitted by Peccci-Quinn strings at T ~ 1 GeV and 
quickly become non-relativistic. The axion contribution to Qq is roughly 

na-ilO~*eY)/ma (71) 

Accelerator experiments plus the requirement that axion emission should not have drained too much 
energy from supernova 1987A require that nia < 10"'^ eV. Since fla < 1, we conclude that if the 
axion exists at all its mass is roughly of order 10"'^ to 10~^ eV, and it gives very roughly of order 
1! One might detect solar axions through their tiny electromagnetic interaction, by exploiting the 
fact that their wavelength (10^*eV/ma) x .197 cm is macroscopic. 

In the context of supersymmetry the axion must have two supersymmetric partners, the axino 
(spin 1/2) and the saxino (spin 0), with rather well defined properties. The former might be a 
dark matter candidate and the latter might dilute the entropy by its decay, reducing predicted 
cosmological mass densities like Eqs. ( [ts]) and ([ti] ) by a factor of up to 10^** [ pT| . 

These dark matter candidates are relatively uncontroversial, involving rather minimal extensions 
of the Standard Model. More ambitious extensions can produce dark matter with practically any 
desired properties; thus, since MDM (mixed, ie hot plus cold, dark matter) has proved so successful 
at fitting the observations, no fewer than three separate proposals have been made for naturally 
generating it in the early universe [^, |l3|, p^ . 

A completely different possibility is that the non-baryonic dark matter consists of black holes 
which form before nucleosynthesis. In that case one has no need of non-baryonic dark matter 
particles. A black hole will form if there is an overdense region with density contrast Sp/p ^ 1, 
whose size is of order the Hubble distance. All of the energy in the overdense regions then ends up 
in the black hole, whose mass is therefore equal to Pr/pm times the mass of the matter in a Hubble 
volume. 



M ^l{)'\prlpm) 



lMeV\^ /lMeV\^ 



T \ T 



10' (72) 



As we shall see later, the density contrast at horizon entry is only of order 10^^ on the scales 
M > lO^^M© explored by large scale structure and the cmb anisotropy, and simple inflation models 
predict that this remains true down to very small scales. The simple models could be wrong however. 
One might also form smaller black holes on scales much less than the horizon size at a phase 



transition. Light enough black holes could be evaporating during nucleosynthesis, or at the present, 
with dramatic effects. 

Finally, I note that although WIMP or axionic dark matter is normally assumed to be gaseous, 
it might be bound into macroscopic objects. To achieve this during radiation domination one needs 
an isocurvature matter density perturbation (one not shared by the radiation) with most of its 
power on small scales. Such a perturbation might be generated during inflation, like the usually 
considered adiabatic density perturbation, but it could also be generated at a phase transition by, 
say, a collapsing bubble. In that case the maximum mass of the objects formed is that within a 
Hubble volume, 

M - .1(1 McY/T^Mq (73) 

This mechanism has been considered in the context of axions, where the phase transition is at 
T - 1 GeV giving objects with M - 10-^° Mq ||. 

Topological defects 

In the early universe, the Higgs fields (and any other scalar fields present) can find themselves in a 
false vacuum, corresponding to a local minimum of their potential. The false vacuum is typically 
stabilized by finite temperature effects, but an interaction with the inflaton field can also to the job 
| p6| , 1^, |l^, ^ |l^. When the temperature, or the inflaton field, falls below some critical value, the 
false vacuum is destabilized and there is a phase transition to the true vacuum. Within the Standard 
Model the electroweak transition is the only example of this phenomenon, but further transitions 
could occur, at much higher energy scales, in extensions of it. 

Such a phase transition can leave behind it topological defects, consisting of regions of space 
where the field is still trapped in the false vacuum. Whether this happens or not depends on the 
nature of the symmetry breaking. No topological defects form at the electroweak transition, but 
they might well form in transitions occurring in extensions of the Standard Model. A concrete 
example is provided by Pcccei-Quinn symmetry breaking, which produces global strings which play 
an important role in producing axions before they disappear at T 1 GeV. Another example is the 
breaking of a GUT which inevitably produces magnetic monopoles. They are so abundant as to be 
cosmologically forbidden unless one either finds a way of annihilating them, or invokes inflation. 

3 The Evolution of the Density Perturbation 

Cosmological perturbation theory develops linear equations for perturbations away from homogene- 
ity and isotropy. Using it one can follow their growth on a given scale, until they become big enough 
for gravitational collapse. On scales > 100 Mpc, where collapse has yet to occur, cosmological per- 
turbation theory can be used right up to the present epoch. On smaller scales the only sure-fire 
way of performing calculations after perturbation theory fails is to perform numerical simulations, 
though analytic approximations can provide some insight. In these lectures I concentrate on the 
linear regime. 

In the Newtonian regime cosmological perturbation theory has long been recognised to be a 
straightforward application of fluid flow equations, taking into account where necessary particle 
diffusion and free-streaming (coUisionless particle movement). In the relativistic regime, cosmological 
perturbations was first discussed by Lifshitz in 1946. His formalism, which has been the starting 
point for most subsequent work including the influential 'gauge invariant' formalism of Bardeen 
considers the perturbed Robertson- Walker metric. An alternative formalism, which makes 
no mention of the metric perturbation and works instead with relativistic fluid flow equations, was 
initiated by Hawking |^3| in 1966. This approach, which treats the Newtonian and relativistic 
regimes in a unified way, is becoming increasing ly popular ||, || ||, H and is the one that I 
will use in these lectures. 

3.1 Relativistic fluid flow 

We populate the universe with comoving observers, who define physical quantities in their own 
region. By definition, the momentum density is zero with respect to a comoving observer; such an 
observer moves with the energy flow. 

A crucial concept is the velocity gradient Uij. It is defined by each comoving observer, using 
locally inertial coordinates in which he is instantaneously at rest, as the gradient of the velocity 



of nearby comoving worldlines, 



p{x,t) = p{t) + Sp{x,t) 
p(x, t) = p{t) + 5p{x, t) 
iJ(x, t) = H{t) + 5iJ(x, t) 



(76) 
(77) 
(78) 



Ui-j = dju^ (74) 

The velocity gradient can be uniquely decomposed into an antisymmetric vorticity ujij , a symmetric 
traceless shear aij, and a locally defined Hubble parameter iJ, 

Uij = HSi-j + cr.y + LUij . (75) 

In the limit of homogeneity and isotropy, Uij = uJij = 0. One can show from angular momentum 
conservation that decays like (p + p)a~^ , so it is presumably negligible. We will see how to 
calculate CTy in Section 3.3. 

Just as for a homogeneous isotropic universe, it is useful to consider 'comoving hypersurfaces', 
defined as those orthogonal the fluid flow worldlines.^ On a given hypersurface, each quantity p, p 
and H can be split into an average plus a perturbation, 



Here t is the time coordinate labelling the hypersurfaces, and x = {x^, x'^, x^) are space coordinates. 
We would like to choose the space coordinates to be comoving coordinates, related to Cartesian 
coordinates by r* = ax"^, with a the average scale factor given by a/ a = H (we shall not have 
occasion to define a perturbed scale factor). This cannot be done exactly, because the expansion is 
not isotropic and the comoving hypersurfaces are not flat. However the departure from flatness is 
of first order, and can therefore be ignored when describing perturbations which are themselves of 
first order. In other words all perturbations 'live ' in flat space. 

Independent scales 

Each perturbation / can be written as a Fourier series, defined in a comoving box much bigger than 
the observable universe 

/(x,t) = ^/k(i)e^'^'' (79) 

k 

The beauty of this expansion is that each Fourier mode propagates independently, as long as the 
cosmological perturbation theory that we are developing here is valid. The inverse wavenumber a/k 
is said to define a scale, which is specified by giving its present value k^^. 

Now consider a small density enhancement in the early universe, which is destined to become, 
say, a galaxy. If its size is of order r = xa, it is typically made out of Fourier components with 
k ~ x~^. As long as it corresponds to a small density contrast, Sp/p ^ 1, it will expand with the 
universe so that its comoving size x remains constant. When its density contrast becomes of order 
1 it will collapse and then its physical size will remain more or less constant. In both cases, though, 
the mass of the enhancement remains fixed. It is therefore useful to associate with each scale k the 
mass of matter enclosed within a sphere of comoving radius x = k~^ (taking the universe to be 
practically homogeneous, corresponding to the early universe). This mass is 

M{x) = 1.16 X 10^^h'^{x/l MpcfMQ (80) 

One expects perturbations with comoving wavenumber k to be relevant for the formation of struc- 
tures with mass M{x), where x = k~^. 

Horizon entry 

The ratio of a given comoving scale a/k to the Hubble distance is equal to aH/k — a/k, which 
decreases with time. At the epoch when this ratio falls through 1, the scale is said to enter the 
horizon. 

Well after horizon entry, the scale is small compared with the Hubble distance, which means 
that ordinary physical effects like diffusion, free-streaming (coUisionless particle movement) and the 
propagation of sound waves can operate, with the expansion of the universe playing only a minor 
role. Well before horizon entry, the scale is much bigger than the Hubble distance, which means 

*As noted later this definition has to be modified if uja does not vanish. 



that causal processes of this kind cannot operate (at least during a Hubble time, and in practice not 
at all). Instead, as we shall see, each part of the universe evolves independently. 
The scale entering the horizon at a given epoch is given by 

k-' = [anr = ^H,-^ (81) 

Except around matter-radiation equality at z ^ 10** one has 

aH oc radiation domination (82) 

aH (X a~^/^ matter domination (83) 

Thus a crude estimate is that the scale entering the horizon at z < 10"* is ^ z^^^'^Hq^, making 
the scale entering the horizon at matter-radiation equality k^^ ^ 10^^ Hq^, and and that the scale 
entering the horizon at z > 10* is fc^^(z) ~ 10^z~^ Hq^ . An accurate calculation shows that k^^ = 
40h~^ Mpc, and that the scale entering the horizon at photon decoupling is k^^^ — 90h~^ Mpc. We 
shall see that the first scale is crucial for structure formation, and the second for the cmb anisotropy. 
The smallest scale directly relevant for structure formation is presumably the one corresponding to a 
dwarf galaxy, which has mass M ~ IO^AIq and from Eq. ( |80| ) corresponds to k^^ ^ .01 Mpc, which 
enters the horizon when T ^ lOkeV. 



The evolution of the density perturbation 

Now I derive differential equations for the perturbations. In doing so one has to remember that 
the comoving worldlines are not in general geodesies, because of the pressure gradient. As a result, 
the proper time interval dr between a pair of comoving hypersurfaces is position dependent. Its 
average may be identified with the coordinate time interval dt, and one can show (using essentially 
the Lorentz transformation between nearby observers) that its variation with position is given by 

^=fl_J^) (84) 
dt \ p + pj ^ ^ 

Along each worldline the rate of change of p with respect to proper time r is given by energy 
conservation and has the same form Eq. (^) as in the unperturbed case, 

^^^3H{p + p) (85) 
dr 

The rate of change of H is given by the Einstein field equation, and to first order receives just one 
extra term in the presence of perturbations, coming from the pressure gradient [ p6[ , 

dH , 4^G IW^Sp 

-7- = -H ^(/3 + 3p)--— — (86) 

dr 6 6 p + p 

This equation is called the Raychaudhuri equation. The operator is the Laplacian on a comoving 
hypersurface, given in terms of comoving coordinates by 

^a-H'=^^ (87) 

Perturbing H , p and p to first order and using Eq. (|^) gives the following equations for the 
Fourier components 

((JpkT = -3{p + p)SHk-3H5pk 

iSH,)- = -2HSH,-^-^Sp. + U'-y^ (89) 

3 3 \a J p + p 

Eliminating (5i?k with Eq. ( p8| ) gives a second order differential equation for pk- It is convenient to 
use the density contrast 6 = in. terms of which the equation is 



H-'Sk+[2-3{2w-ci)]H-'Sk--il~6ci + 8w-3w')S^^-l — ] (90) 



I have used the notation w ^ p/ p and = pj p. 

In these equations for the perturbations the distinction between p and p is not meaningful 
(similarly for p and H), because we are working only to first order. For the same reason the 
distinction between r and t is not meaningful. In order to have linear equations one will therefore 
obviously take p to mean p and and similarly for p and H , and will take the dot to mean differentiation 
with respect to t. Note that = pj p is the speed of sound, because p and p vary adiabatically in a 
homogeneous isotropic universe (no heat flow). 



The case of zero pressure gradient 

The right hand side of Eq. (pO|), which involves the pressure gradient, is negligible after matter 
domination because the pressure is negligible.^ It is also negligible well before horizon entry even 
during radiation domination, because the gradient is so small.^ When it is neglig ible, Eq. (^) can 
be reduced to a first order equation, which has a very simple interpretation p^ , 

The solution and its interpretation hinge on the introduction of a quantity defined locally 
through the Friedmann equation Eq. (|^). Just as in the homogeneous, isotropic case, general rel- 
ativity shows that K/a? is a measure of the curvature of the comoving hypersurfaces. In the 
homogeneous, isotropic case K is time independent by virtue of the energy conservation Eq. ( ^ ) 
and the gravitational deceleration equation Eq. (^). We have seen that in general the latter is 
modified to become the Raychaudhuri equation Eq. (p6|), and as a result K is not in general time 
independent. But when the pressure gradient is negligible they are both unmodified, so that K is 
time independent. We can say that when the pressure gradient is negligible, each region of space 
evolves like a separate Friedmann universe. 

On a comoving hypersurface K can be split into and average K plus a perturbation SK, but the 
average can be set equal to zero because i7 ~ 1. Perturbing the Friedmann equation therefore gives, 
to first order, 

2HSH, = '-^Sp, - ^ (91) 
3 

When SKy^ is time independent, Eqs. (|8^) and ( |9l| ) give a first order differential equation for the 
density contrast, 

2H-^ d 



5 + 3w dt 



aH^ ^ 



5 



oHV „ 2 + 2w 



where w = p/ p and I have introduced the useful quantity 



Remembering that SK/a^ is the curvature perturbation and that it has units of (length)"^, we 
see that TZ\^ — (9/e)((5/C/H^)(H^/|j'=) essentially measures the curvature perturbation i n un its of 
the relevant scale a/k. (The factor 3/2 is chosen so as to give the simple relation Eq. ( |l99|) with 
the inflaton field perturbation.) As we shall verify shortly, another interpretation of TZ is that it is 
essentially the Newtonian gravitational potential caused by Sp. 

During any era when w is constant, Eq. (^2|) has the solution (dropping a decaying mode) 

aHV ^ 2 + 2w 



In the radiation dominated era before horizon entry this becomes 



fc y 9 

and in the matter dominated era it becomes 

oHV ^ 2 



^) (5k = ^7^k (initial) (95) 



— j <5k = -7^k (final) (96) 



^Except for the baryons on scales below the Jeans scale, and we are presently assuming that the dark matter is 
cold, otherwise it is modified by free-streaming. Both of these points will be addressed shortly. 

^"Provided that p/p is not extremely large, which is ensured by the adiabatic initial condition defined shortly. 



As the labels imply I am regarding the value of 7?.k during the first era as an 'initial condition', 
which determines its value during the 'final' matter dominated era. 

For future reference note that during matter domination, H (x t^^ oc a^^^^ and 

oc a (matter domination) (97) 
3.2 The transfer function 

For scales entering the horizon well after matter domination {k^^ ^ k^^ = 40/i"^ Mpc) the initial 

and final eras overlap so that 7?.k(initial) = 7?.k(final). Otherwise there is a transfer function T{k), 
which may be defined by 

Tek (final) = T(fc)7^k (initial) (98) 
An equivalent, and more usual, definition is 

a"Mk(final) = ylT(fc)<5k(initial) (99) 

where the (time dependent) right hand side is evaluated at an arbitrarily chosen time during the 
initial era, and the constant A is chosen so that T becomes equal to 1 on large scales. 

The adiabatic initial condition 

In order to calculate the transfer function one needs an initial condition, specifying the density 
contrast of each species of matter and radiation before horizon entry and before matter domination. 
The matter consists of baryons and one or more species of non-baryonic dark matter. The radiation 
consists of photons and massless neutrinos. 

The most natural condition, is the adiabatic condition, which is that the density of each species 
depends only on the total energy density, or equivalently on the temperature. In other words, each 
density contrast vanishes on a hypersurface of constant energy density. Going along any comoving 
worldline, the density pr of any species of radiation varies like and the density pm of any species 
of matter varies like a~^, so that 

3 

Spm/P^n = -Spr/Pr (100) 

Thus, the adiabatic initial condition implies that each species of radiation has a common density 
contrast and so has each species of radiation, and that they are related by a factor 3/4. 

The most general initial condition is a superposition of the adiabatic one and some isocurvature 
initial condition, specifying a set of initial density perturbations 5pi which add up to zero. An 
isocurvature initial condition is not very natural, especially in the context of inflation. Furthermore, 
it turns out that a purely isocurvature intial condition (with a reasonably flat, Gaussian spectrum) 
cannot lead to a viable structure formation scenario. The adiabatic condition is assumed from now 



Cold and hot dark matter 

The other ingredient needed to calculate the transfer function is the nature of the dark matter. It is 
practically always taken to have negligible interaction with anything else (to keep it dark). Neutrino 
dark matter is called hot dark matter (HDM) because it remains relativistic (in other words it is 
radiation rather than matter) until a very later epoch, which is after structure forming scales have 
entered the horizon for pure HDM. Cold dark matter (CDM) by contrast is nonrclativistic at all 
epochs of interest, and warm dark matter is in between. Only CDM and HDM have as yet been 
studied in the light of modern data. Pure HDM is ruled out (at least if we do not form structure with 
seeds such as cosmic strings, or with a density perturbation having so much small scale power that 
very small objects form first). Pure CDM performs rather impressively, going a long way towards 
explaining several different types of data, but it cannot actually fit them all simultaneously within 
their accepted uncertainties. For this reason people have considered mixed dark matter (MDM), 
roughly 70% cold and 30% hot, which at the time of writing seems to do rather well pG, ra, M. 



Figure 2: The transfer function T{k), taken from Pogosyan and Starobinsky | |3C| |. From top down, 
the curves correspond to Vli, = 0, 0.05, ...0.95. The transfer function is evaluated at the present 
epoch, but in the hmit of pure CDM (top curve) it becomes time independent. 



The CDM transfer function 

The transfer functions for pure CDM and for MDM are shown in Figure 2. Although several 
different effects have to be taken into account to calculate them accurately, their general form is 
easy to describe. Let us deal first with the pure CDM transfer function. 

On a given scale essentially nothing happens before horizon entry because there is no time for 
particles to move on that scale. As a result, the growth 5 cx (aH)"^ of the density contrast, 
that we established earlier as a consequence of the 'separate Friedmann universe' evolution of each 
comoving region of space, is shared by each species. After horizon entry, the density contrast of 
massless neutrinos rapidly decays away, because the neutrinos free-stream (that is, travel freely) 
out of any over dense region. The baryons and photons are thermally coupled until the epoch of 
photon decoupling, and their density contrast oscillates as a sound wave, which decays more or less 
rapidly (depending on the scale) as the photons diffuse away from overdense regions. The density of 
cold dark matter grows, but until matter domination the amount of growth is rather small because 
the matter does not have much self-gravity. For example, taking the radiation to be uniform and 
ignoring the baryons, one finds by following the steps leading to Eq. ( |9C| ) that the matter density 
contrast satisfies 



Using Pm/ Pr = a/ acq, one finds that the solution of this equation is proportional to 2 + ia/ acq (plus 
a decaying mode), which indeed hardly grows before matter domination at a = acq. Since it would 
have grown like (aH)~^ had horizon entry not occurred, the total amount of growth missed by a 
given scale is roughly 



5™ + 2H5, 



m 



(101) 




(102) 



The CDM transfer function is therefore roughly 



T{k) ^1 
T{k) c 



(fc ^ > kcq - 
(fc-1 < fc-1). 



40/1-1 Mpc) 



(103) 
(104) 



This indeed gives an extremely crude representation of the correct CDM transfer function. To 
calculate it accurately one needs to follow in detail the oscillation and/or decay of each component.^ 

The above discussion applies only to the CDM density perturbation. We need to know also what 
happens to the baryon density perturbation, since baryons after all are what we observe. Unlike the 
CDM density contrast, the baryon density contrast is small at the photon decoupling epoch because 
it has been decaying since horizon entry. Afterwards, the baryons are unaffected by the photons, 
but two competing forces act on them. First there is gravity, which tends to make them fall into the 
potential wells caused by the CDM density contrast, and second there is their own pressure gradient 
which tends to keep them out of the wells. 

To rigorously see which effect wins one should generalize Eq. ( pO[ ) to treat the CDM and the 
baryons as a pair of uncoupled fluids ^ . In practice the following order of magnitude estimate 
is enough. Ignoring the pressure, the time taken for the baryons to fall into a well isj^ of order 
(Gp)~^/^. The time taken for the pressure to adjust itself to prevent the collapse is of order A/cs 
where A = 2i:/k is the wavelength and Cs is the speed of sound. Collapse occurs if A/cs < (Gp)"^/^ 
because the pressure cannot act quickly enough. Inserting by convention a factor 2, one concludes 
that collapse occurs on scales in excess of /cj = {^.TrGp/vlY^^ . This is called the Jeans scale and the 
corresponding mass, given by Eq. (^), is called the Jeans mass. 

Just after decoupling the speed of sound is given by = ST/Stoat where T ~ 1100 x 2.74K. 
The corresponding Jeans mass is of order IQ^Mq. As yet structure formation on such small scales 
is not well understood because the necessary numerical simulations of the collapse are too difficult, 
but it may be significant that there are no galaxies with mass less than the Jeans mass. 



The MDM transfer function 

The hot (neutrino) component becomes non-relativistic only at the epoch 

1 + z„,. = (1.7 X 10-^)/i2rj^ (105) 

On scales entering the horizon before this epoch, the neutrino dark matter will free-stream out of any 
density enhancement so that its density contrast will decay. Afterwards the HDM density contrast 
is free to grow to match that of the CDM, on scales in excess of an effective Jeans length 

fcj^ = .ll(l + z)i/2f7^Mpc (106) 

By the present time, the CDM and HDM have a common density contrast on scales > .1 Mpc, 
which are the only ones of interest. Note, though, that the MDM continues to evolve right up to 
the present. 



3.3 The peculiar velocity field 

After matter domination it is very useful to introduce the concept of a peculiar velocity field. 
Traditionally this concept is formulated in the context of the Newtonian treatment of perturbations, 
but it has been recognised recently that the same concept is valid, and useful, also in the general 
relativistic treatment ||^, |3l| . 

At a given epoch the Newtonian treatment is valid on scales much less than the Hubble distance, 
to the extent that gravity waves can be ignored. There is a well defined fluid velocity u, and choosing 
the reference frame so that u vanishes at the origin the PV field V is defined^ by 

u(r) - u(0) = Hr + v(r) - v(0) , (107) 

An equivalent statement in terms of the velocity gradient Eq. (u5^ is 



5uij — diVj . (108) 

^^A surprising V good approximation ]32| is to ignore the decay, treating the radiation and matter as uncoupled 



perfect fluids 
^^If a partic 
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but I do not know of any simple reason why this is so. 
Tails from rest towards a point mass M , its velocity at distance r is given by mv'^ = 2GM/r so it falls 
a signiflcant distance in a time t ~ r/i) ~ (GAf/r»)-i/2. We are replacing the point mass by a perturbation with size 
r and density p ~ M/r^. 

^^Up to a constant, which can be chosen so that the average of v vanishes. 



where di = a~^d/dx'^ . Like any vector field, v can be written v = v'^ + v""", where the transverse 
part v""" satisfies divj = and the longitudinal part is of the form v'^ = VV'ti- Eq. defines the 
local Hubble parameter shear aij and vorticity Wy , this last being given by 

^^3^\{^^V] -d,vj) (109) 

Angular momentum conservation gives uoij — s- a^^, so v"^ decays like and may be dropped 
(remember that di = a~^d/dx^). 

Taking v to be purely longitudinal, it is determined by the density perturbation in the following 
way First, take the trace of Eq. (|l|) to learn that V.v = iSH. From Eqs. (H), (|l|), (^ and (|9|), 
it follows that 

V.v = -(47rG(5p)i (110) 

The solution of this equation is 

^^=-tV'^p (111) 

or 

v,{^,t)^-{tla)^^ (112) 
V'(x) = [ ^-f^^dPx' (113) 



where 



The factor converts coordinate distances into physical distances. Since it is related to the 
density perturbation by the Newtonian expression, ip is called the peculiar gravitational potential. 
It is independent of t because, from Eq. (p7|), 5p a? . 

From Eq. ( p6|) we see that the peculiar gravitational potential is related to the spatial curvature 
perturbation by 

ij = -|7^(final) (114) 



From Eqs. (Ill) and (113) the Fourier components of v, ij: and 5 are related by 



3 f aH"- ^ 



- (116) 

The extension of these results to the relativistic regime turns out to be very simple The 
vorticity decays (p + p)a~^ and is again presumably negligible. In that case one can show that the 
velocity gradient perturbation is of the form 

Suij = + i/ijj (117) 

The extra term, which is absent in Newtonian physics, represents the effect of gravitational waves. 
The success of Newtonian theory on scales < 100 Mpc implies that it is negligible there, but as 
we shall see it could be significant on larger scales, and contribute to the cmb anisotropy. Even if 



= 0. 



present though, it does not spoil Eqs. (110)-(116) because it is transverse, dihij = 0, and traceless. 



One can avoid dropping the vorticity by proceeding as follows [pll. First, a transverse velocity v 



may be defined by Eq. (109). The worldlines with velocity — relative to the comoving worldlines 
have no vorticity, and comoving hypersurfaces are defined to be orthogonal to them (there are no 
hypersurfaces orthogonal to worldlines with nonzero vorticity). The velocity gradient 6uij receives 

an extra contribution \ {divj — djvj) + ^{diwj + djwj) where w'^ = 1 + Q {l + ^ (t)^ ^'^ ^^"^ 
scales well inside the horizon w is negligible and the Newtonian result is recovered. 



3.4 The spectrum of the density perturbation 



In order to discuss the perturbations in a given region of the universe around us, one has to perform 
the Fourier expansion Eq. (|7^) in a box much bigger than this region. If the box is a cube with sides 
of length L, the possible values of k form a cubic lattice in k space with spacing 2tt/L. 

When discussing an isolated system, which is the usual case in physics, one can take the limit 
L — > (X) in a straightforward way, the coefhcients /k tending to a constant limit which is a smooth 
function of k. But cosmological perturbations do not fall off at large distances, and their Fourier 
coefRcients are not smooth functions of k. They are the spatial analogue of a signal extending over 
an indefinite period of time, as opposed to an isolated pulse. 

Although the coefRcients /k are not smooth functions, it is reasonable to suppose that |/kp is 
smoothly varying when smeared over a region d^fc of fc space, which is large enough to contain many 
lattice points. I shall denote this average by (|/kP)- It depends only on A: = |k|, and up to a A: 
dependent factor it is called the spectrum of /, because of the analogy with a signal. A convenient 
choice of the factor is to define the spectrum as 

Vf ^ An{\h\') (118) 

The normalisation is chosen to give a simple formula for the dispersion (root mean square) of /, 
which I shall denote by (T/. From the Fourier expansion one has a'j = X) l/kl' ^^'^ since the possible 
values of k form a cubic lattice with spacing 2tt / L the transition from sum to integral is 



(2-K 



^ — > 47r / k'^dk (119) 



The dispersion <J/ is therefore given by 

aj^{fi^))^J%fik)^ (120) 

with the brackets now denoting the average over position x. 

For the density perturbation / = J it is useful to define the correlation function ^(r) by 

f(r) = (/(r + x)/(r)) = J^^ Vf{kf^^f (121) 

The analogous quantity is useful for other perturbations like the peculiar velocity components, 
though it is not then called the correlation function. For r = it clearly reduces to cr^ . 

If the phases of the Fourier coefficients are random, / is said to be Gaussian, and then all of its 
stochastic properties are determined by its spectrum. In particular the probability distribution of /, 
evaluated at randomly chosen points, has a Gaussian profile. As long as the perturbations are evolv- 
ing linearly this Gaussian property is true of perturbations originating as inflationary fluctuations, 
and I take it for granted from now on. 

From Eqs. ( |9q ) and (p8|), the spectrum of the density contrast after matter domination may be 
written ^ 

-Psik) = T\k)5Uk) (122) 



The quantity 5h specifies the initial spectrum. In fact, from Eq. ( |95D it is related to the spectrum 
of the initial curvature perturbation TZ by 

5Uk) = ^Vn{k) (123) 

The subscript H has been used because 51^ is exactly equal to the value of Vs on horizon entry on 
scales fc^^ ^ fcjjq^, and approximately equal to it on smaller scales. As we shall see later, this means 
that Snik) is roughly equal to the rms value of 5 at horizon entry, for a density enhancement with 
comoving size of order k^^. 

The standard assumption is that 5^ is independent of k. A more general possibility is to consider 
a spectrum 

5]j cx (124) 



where the exponent n is cahed the spectral index. (The definition of the index as n — 1 instead of n 
is a historical accident.) The standard choice of n = 1 was first advocated by Harrison (1970) and 
Zel'dovich (1970) on the ground that it is the only one making the perturbation small on all scales, 
at the epoch of horizon entry. Although this is a powerful argument in favour of a value of n of 
order 1, it does not require that n is actually close to 1 unless one assumes, without justification, 
that the power law dependence holds over a huge number of decades of scale^ 

In this context, infiation makes a striking prediction p3| . It gives the spectral index is given in 
terms of two parameters ei and ?7i , 

n = 1 + 2771 - 6ei (125) 

The parameters depend on the model of inflation, but both are much less than 1 in magnitude and 
ci is positive. As a result most inflationary models predict that n is rather close to 1, and typically 
smaller rather than bigger. The prediction that n is close to 1 seems to be conflrmed by observation, 
which is a striking success for the general idea of inflation. In the very near future, one will be able 
to use the observed value of n as a powerful discriminator between different inflationary models. 



4 The Cosmic Microwave Background Anisotropy 

The first detection last year of the intrinsic anisotropy of the cmb was surely the most important 
advance in observational cosmology for a long time. Exploring scales two orders of magnitude bigger 
than the biggest that can be probed by galaxy surveys, the fact that it agreed within a factor of 
two with extrapolation from these surveys using the pure CDM model with a flat spectrum was 
extremely impressive. Before exploring the significance of this fact in the next section, let us study 
the cmb anisotropy is some detail. 



4.1 The spectrum of the cmb anisotropy 

At a given wavelength, the cmb is anisotropic ie its intensity depends on the direction of observation. 
The anisotropy is usually specified by giving the equivalent variation in the temperature of the 
blackbody distribution. Denoting the direction of observation by a unit vector e, the anisotropy 
may be expanded into multipoles 

^ + £ E ^^rie) (126) 

1=2 m=-l 

The dipole term w.e is well measured, and is the Doppler shift caused by our velocity w relative 
to the rest frame of the cmb (defined as the frame where it has zero momentum density). Unless 
otherwise stated, AT will denote only the non-dipole part from now on. Its mean square over the 
whole sky is 

1=2 m=-l 

The multipoles, and therefore the mean square anisotropy, depend on the position of the observer. 
Averaging over position gives the result for a randomly placed observer, 

1 °° 

= i^E(2^ + l)^' (128) 

position 1 — 2 

where 

^' = {\<n,os-.:on (129) 

In contrast with the quantities Tjf are expected to be smooth functions of For large I 

one can therefore replace the sum by an integral, 






position 



-I poo TJ 

-J^ m^mff (130) 



^"^On cosmologically interesting scales the observed magnitude of 5h is of order 10 ^, so if n = 1 — (5n ~ 1 one has 
to extrapolate over 5/Sn orders of magnitude before Sjf becomes of order 1. 



Figure 3: This figure is reproduced from |34 . The top Figure shows the spectrum of the cmb 
anisotropy for a shghtly tilted standard CDM model with spectral index n — 0.85 and Hb = 0.05. 
As discussed in the text, the case n — 1 is recovered for Z ^ 1 by multiplying the curves by a 
factor — l^'^^. The contribution of an adiabatic density perturbation is the middle line, 

labelled 'S', and the contribution of gravitational waves is the bottom line, labelled 'T'. The light 
dashed line is the density contribution with fi^ = 0.01. For each curve, the quantity plotted is 
1(1 + , normalised to 1 for the quadrupole, 1 = 2. If the gravitational wave contribution to 
the quadrupole is equal to that of the density perturbation, as is roughly the case for power-law 
inflation with n = 0.85, the top curve indicates the total. On the other hand, it could well be 
that the gravitational contribution is negligible. The bottom Figure shows the filters Fi for various 
experiments, which are cited in ||3^ . 

Another nice feature of large I is the relation 6^1/1 between the angular size of a feature in 
the sky (in radians) and the order I of the multipoles that dominate it. (This is analogous to the 
relation r ~ between the linear size of a feature and the wavenumbers that dominate its Fourier 
expansion.) Translating to degrees we have the following relation between I and the angular scale 

GO , ^ , 

To - - (131) 



By analogy with the Eq. (120) for the density contrast, one may call (2Z-I- l)ZS^/47r the spectrum 
of the cmb anisotropy. 

4.2 Predicting the cmb anisotropy 

The predicted spectrum of the cmb anisotropy is shown in Figure 3. I will discuss the prediction 
qualitatively, then give a detailed calculation of the large scale (small I) part which is the most 
interesting. 

Discounting the possibility of early re-ionisation the cmb last scattered at the epoch of photon 
decoupling. Its surface of last scattering therefore lies practically at the particle horizon whose 
comoving distance is x = 2i?Q"^ = GOOO/i""'^ Mpc. At this surface an angle 9 degrees subtends a 
comoving distance 

X 6 

lOO/i-iMpc "1° ^ ' 

This gives, at least roughly, the linear scale explored by observing the cmb on a given angular scale. 



It is convenient (at least pedagogically) to separate the anisotropy into an initial anisotropy, 
specified on a comoving hypersurface shortly after last scattering, and the additional anisotropy 
acquired on the journey towards us, 



Ar(e) /AT{e)\ f AT{e 



T \ T J \ T 

\ / cm \ 



(133) 



The initial anisotropy represents the variation of the intensity of the cmb, as measured by different 
comoving observers, each with their detector pointing in the relevant direction e. 

Let us discuss these contributions, from the largest scales downwards. Scales in excess of a degree 
or so explore linear scales in excess of 100 Mpc, neatly complementing galaxy surveys which explore 
scales < 100 M pc. S uch scales have not entered the horizon at last scattering, so the adiabatic initial 
condition Eq. ( |100| ) still holds at that epoch. The cmb is practically isotropic for each comoving 
observer, but its intensity varies from place to place because pr oc T^. Since last scattering occurs 
during matter domination, the density contrast is practically that of the matter, so the adiabatic 
condition gives 

' ^ ^ ^ - :'^(Xem,tcm) (134) 



V T 3 

However we shall see that this effect is dominated by the additional anisotropy acquired on the 
journey towards us which is called (on these large angular scales) the Sachs- Wolfe effect. 

Scales in the arcminute regime explore distance scales of order 10 to 100 Mpc, the same as the one 
explored by large scale galaxy surveys. As the scale decreases through this regime, the anisotropy at 
first exhibits a peak, which can be interpreted as of the Doppler shift caused by the peculiar velocity 
of the initial comoving hypersurface. Thereafter the anisotropy acquired on the journey towards us 
becomes negligible]^ The initial anisotropy also falls, but it remains significant until the scale falls 
to a few arcseconds. 



On scales less than a few arcminutes, the distance Eq. (132) subtended at the last scattering 
surface is less than the thickness of this surface. As a result the cmb is isotropised by its passage 
through the last scattering surface, and it turns out that no significant anisotropy is acquired by the 
cmb on its journey towards us either. The anisotropy is therefore expected to be very small. 

The cmb transfer function and the cosmic variance 

By evaluating these physical effects, one can calculate the multipoles a™ seen by an observer at any 
given position, in terms of the initial density perturbation (5k (initial).^ Of course we do not know 
(5k(initial), nor are we interested in it. Instead we know, or rather are prepared to hypothesise, its 
stochastic properties, which are that it is Gaussian with a spectrum specified by S'jj{k) cx fc"^^. The 
Gaussianity implies that the real and imaginary part of each multipole a™ has a Gaussian probability 
distribution as a function of the observer's position, with no correlation between them (except for 
the trivial one aj"* = fl;"™). The spectrum Yif of the cmb anisotropy — specifying the widths of the 
Gaussian distributions — is related to the spectrum Snik) of the initial density perturbation by a 
transfer function Ti(k), 

Ef= r T?{k)5l{k)'^ (135) 



Jo 

In accordance with Eq. ( |132| ), T;(fc) is peaked at 

I 100/i-iMpc 

60 " fc-i ^^^^^ 

The transfer function is independent of the nature of the dark matter on large scales as we discuss 
in a moment. On smaller scales the effect of adding in some hot dark matter has not yet been 
investigated in detail, but going from pure CDM to MDM (30% hot dark matter) is not expected 
to have a big effect. The prediction in Figure 3 is for pure CDM. 

This prediction gives the expectation value of the multipole |a™p evaluated at a random 
position, but we can only measure |a™p at our position. The best guess is that |a™p = E^, but one 

^^I here discount the Sunyaev-Zel'dovich effect by which the cmb acquires anisotropy through photon up-scattering 
off dust on its way through a galaxy or cluster, ft occurs on scales of order one arcminute and is easily identified for 
nearby clusters. The cumulative effect of distant clusters is expected to be negligible. 

^^Discounting gravitational waves, which f consider in a moment, and isocurvature density perturbations, cosmic 
strings etc which I do not consider at all. 



can also calculate the expected error, which is the square root of the variance of the guess (the mean 
value of (|a™P — T,f)'^). This variance is called the cosmic variance. For the real or imaginary part of 
a single multipole, the Gaussian distribution implies that the variance is equal to 2 times the mean 
square prediction. The Ith multipole has 21 + 1 independent components, so summing over m makes 
the cosmic variance only 2/(21 + 1) times the mean square prediction. A typical observation typically 
involves a (weighted) sum over several multipoles, which further reduces the cosmic variance so that 
in practice it is not usually a dominant factor. One must always bear it in mind though. 



4.3 The Sachs- Wolfe effect 

Now we calculate the anisotropy on scales in excess of a few degrees ^ . 

Consider a photon passing a succession of comoving observers. Its trajectory is adx/di = — e 
and between nearby observers its Doppler shift is 



d-^ , da . , 

— = eiCjUijclr = — — + eiejOUijdr , (137) 
A a 



where the first term is due to the average expansion, and the second is due to the relative PV of 
the observers. Integrating this expression gives the redshift of radiation received by us, which was 
emitted from a distant comoving source. The unperturbed result is A/Acm = l/ocm, and the first 
order perturbation gives 

AT(e) /AT(e)\ 

-^=[-^) e,e,Su,,{^,t)a{t)dx , (138) 

where Xem — "^H^^ is the coordinate distance of the last scattering surface, and the integration is 
along the photon trajectory 

xit) = / _ = 3 p - - . (139) 
Jt a \ao aj 



Using Eqs. ( |108D , ( p^39D and (111) and integrating by parts one finds 

e- [v(0,io)-v(xe,„,tom)] . (140) 

Here Xem = a^emG is the point of origin of the cmb coming from direction e. 

A better expression follows |^ if one uses the divergence theorem and Eq. (13£) to project out 
the dipole part of '0(xom)/3. One finds that it is equal to (v(tom) — v(io)) where () denotes the 
average within the last scattering surface x = Xom- Defining v = v — (v) this gives 

Ar(e) /Ar(e)\ ^1.,, ^ 

J. = { T ) + 3 ['^(^<='°)]i>i + 

e- [v(0,io)-v(xe„„tcm)] • (141) 

On angular scales ^ 1°, corresponding to linear scales at last scattering which are outside the 



horizon, Eqs. ( |115| ) and (116) show that 

IV'kl » I'^kl » |vk| (142) 
On these scales we can therefore drop the term v(xc,n,iom)j as well as the the initial anisotropy 



Eq. ( |134| ). This leaves 

^^ = i[V.(xe„,)],>i+e.vo. (143) 

The last term in this expression is the dipole and it defines the rest frame of the cmb by giving 
our velocity vo = vo — (vo) relative to that frame. Since vo is our peculiar velocity, the peculiar 
velocity of the cmb rest frame is (vq), the average peculiar velocity (bulk fiow) of the region inside 
the surface of last scattering (at roughly the Hubble distance from us). As we shall discuss in Section 
5 the bulk flow decreases as the scale increases, and is quite negligible on the Hubble scale, so in 
practice we can just say that the rest frame of the cmb has zero peculiar velocity. This indeed is 
what is usually assumed. 



Inserting the Fourier expansion of ?/; and projecting out a multipolc leads, for an observer at the 
origin of coordinates, to 



(144) 



where ilk is the sohd angle in k space. For a randomly placed observer the coefficients (5k have 
random phases (except for the trivial correlation implied by the reality condition = (5-k), and 
this implies that the multipoles a™ have a Gaussian distribution. The variance of the distribution 
is (Peebles 1982b) 

Sf^^y^ -jf{2k/aH) 5l{k) (145) 
where ji is the spherical Bessel function. With Sjj{k) cx fc"^^ this becomes 



r((3-n)/2)r(Z + (n-l)/2) 



r((4-n)/2)r(Z + (5-n)/2) 



l{l + l) 



(146) 



The square bracket is equal t o 1 f or n — \. For I ^ 1 and / ^ |ri| it can be replaced by 1, if Snik) 
is evaluated on the scale Eq. (136) which dominates the integral. 



4.4 The contribution of gravitational waves 

We have so far ignored gravitational waves. They correspond to a metric perturbation hij which is 
traceless, S^^hij = 0, and transverse, dihij — 0. This means that each Fourier component is of the 
form hij — h^efj + hxe^j, where in a coordinate system where k points along the z-axis the nonzero 

e+ = 1 and e^y = e^^ = 1. The 



components of the polarisation tensors are defined by e^^ = 
spectrum Vg of the gravitational wave amplitude may be defined by summing Eq 
components, 

^^^KS) ((l^+(k)P> + (|/^x(k)p)) 



3^ over all four 



(147) 



Each Fourier component satisfies the massless wave equation, which in comoving coordinates is 



[k/aYhi 







(148) 



Well before horizon entry it has constant initial value. For scales entering the horizon after matter 
domination its subsequent evolution is 



hij (initial) 



(149) 



where x = 2k/ (aH). Well after horizon entry one has redshifting radiation. 

As noted in Eq. (117), the contribution of gravitational waves to the velocity gradient is 6uij — 



hij/ 2. By substituting this expression into Eq. ( |138| ) one can calculate the cmb multipoles in terms 
of the initial amplitude, and hence calculate cmb spectrum terms of the initial gravitational 

wave spectrum Vg{k). Each gravitational wave gives its dominant contribution as it enters the 
horizon, since its amplitude is practically constant before that and redshifts away afterwards. As 
a result the gravitational wave contribution to the cmb anisotropy cuts off below the angular scale 
~ 1°, corresponding to the smallest linear scale which enters the horizon after decoupling k~^ — 
90/i^^ Mpc. The corresponding multipole cutoff is at Z ~ k^^^Ho/2 ~ 70. 

On scales well in excess of a degree {I < 10), the gravitational wave contribution is given in the 
case of a flat spectrum (Vg(k) independent of k) by Starobinsky pst . 



dec 
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(150) 



If one ignored the cutoff due to the redshift, the coefficient C; would become 1 in the limit ^ 3> 1. 
Starobinsky gives the values C2 = 1.118, C3 = 0.878 and C4 ~ 0.819. Numerical calculation, 
including the effect of the cutoff, shows that a value C; close to 1 is indeed achieved for I ^ 10, 
before the cutoff takes effect (Figure 3). 



For / ^ 1 the above result is good also iiPg{k) has moderate scale-dependence, provided that it 



is evaluated at the scale Eq. (136) which dominates the ^th multipole. 

Within a given inflationary model one can calculate Vg{k). Defining the spectrum ng of the 
gravitational waves by Vg oc fc"9 , one flnds 

Ug = -2ei (151) 



where ei is the same small positive parameter appearing in Eq. (12f;) for n. Setting the coefficient 
Ci in Eq. (|l50| ) equal to 1, the ratio of the gravitational and density contributions is also given in 
terms of this parameter, by [B3l 



S^(grav) ^ 

^=S2(density)-^^'^ ^^^^^ 



In some models of inflation, ei is very small, corresponding to a very small gravitational wave 
contribution with a very flat spectrum. If ei is signiflcant, the potentially is typically of the form 
with a at least equal to 2 and often much bigger (I include the behaviour e"*"^ in this class since 
it corresponds to the limit a — > oo). In that case one has — ?ig = 1 — rt > and 

i?~6(l-n) (153) 

Thus one's expectation from inflation is that either gravitational waves are negligible or that 



their relative magnitude is related to the spectral index by Eq. (153) (on scales in excess of a few 
degrees). It should be noted that this expectation comes from particle physics, being a consequence 
of the kind of potentials that typically arise. It is not generic to the idea of inflation per se, which 
in fact provides no relation whatever between R and n since the two parameters ei and 771 can be 
chosen independently. 

4.5 Observing the cmb anisotropy 



By analogy with Eq. ( 121 ), one can define a temperature correlation function, 



Here e and e' specify the directions in which the anisotropy is observed, and the average goes over 
directions separated by an angle a. It is given in terms of the multipoles by 

00 

C{a) ^^QjPi (cos a) (155) 

where 



^ i E (156) 



-I 



Since one is not interested in very small scale fluctuations, it is convenient to smear AT(e)/T 
over a patch of sky with angular size 9f, to obtain a 'filtered' anisotropy AT{0 f,e)/T. This cuts 
off the multipole expansion above I ^ ^J^y to be precise, with a Gaussian profile for the smearing 
function, its effect is 

ar^exph(2Z + l)6»//2]ar (157) 
Associated with the smeared quantity is a correlation function 

00 

C{9f,a) = ^exp [-{{21 + 1)6^/2)^] QfPi{cosa) (158) 

A given experimental setup typically measures something which can be directly related to a 
suitably smeared C. The simplest one is a single beam whose resolution can be represented by a 
Gaussian profile. Averaged over the sky the anisotropy measured by such a beam is 



aUOf) 



T 



C{0f,O) (159) 



For more complicated setups involving two or three beam switching, still with Gaussian profiles, the 
measured anisotropy is 



mOf^e) AT{9f,e') 



T 



T 



(160) 



and 



AT(gy,e) 
T 
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^C(0/,O)~2C(0/,a) + ic(0/,2a) 



(161) 



In the second expression, e' and e" lie on opposite sides of e, aligned on a great circle and each at 
an angular distance a. In a typical setup the 'beam throw' a is of the same order of magnitude as 
the 'antenna resolution' 9f. 



If we denote the left hand side of Eq. (15£), (160) or (161) generically by (AT/T)^, the prediction 
in terms of multipoles may be written 



oo 



1=2 



(162) 



where Fi is a filter function. For the single beam expression Eq. (15£) the filter keeps all multipoles 



with I <9j , but for the two- and three beam expressions it keeps only multipoles with I ^ 9 

For a realistic experimental setup these Gaussian expressions require modification. Accurate 
filter functions for some currently mounted observations are shown in Figure 3, which is reproduced 
from lal. 



Present observations 

Until recently the only published observation of the non-dipole cmb anisotropy was that of the 
COBE group |36|. They give the quadrupole and some information about higher multipoles, but 
because of cosmic variance their most useful statistic is the rms anisotropy on the scale 10", the 
smallest that they explore. 

aT{9f) = (1.10 ± .18) X 10"^ (163) 

Assuming a Gaussian beam profile one can use this expression to normalize the spectrum Snik) of 
the density perturbation, given the spectral index n and the relative contribution R of gravitational 
waves. With n — 1 and i? = it gives Sh — (1.7 ± .3) x 10~^. Subsequent analysis |3^ using the 
correct beam profile and allowing for incomplete sky coverage raises this to by about one standard 
deviation to 

5h ^ (2.0 ±. 3) xlQ-^ (164) 

Since the ten degree measurement probes scales of order the Hubble distance k^^ ~ lO^Mpc, 
this provides a reasonable estimate of Sh on that scale even if n ^ 1 is somewhat different from 1, 
provided that gravitational waves are negligible. Most inflationary models predict either that they 
will be negligible, or t hat they will be given by Eq. ( 153 ) (with n < 1). In the latter case the 
normalisation Eq. (164) becomes 

Sh = (2.0 ± .3) X 10-^[1 + 6(1 - n)]"^/^ (165) 

At the time of writing additional measurements on smaller scales are being reported. They have 
bigger uncertainties than the COBE measurement, and do not add substantially to what we learn 
by combining this measurement with the results of galaxy surveys in the manner discussed in the 
next section. 



5 The CDM and MDM Models of Structure Formation 

A model of large scale structure should simultaneously account for at least half a dozen completely 
different types of data, exploring a wide range of scales. Pure CDM with a flat spectrum (n = 1) 
performs rather impressively, but cannot actually fit all data within the accepted uncertainties, and 
the same remains true if one allows n to vary. Adding about 30% of hot dark matter does work, as 
far as theory and observation have been compared 



5.1 The filtered density contrast 



At the present epoch the universe is highly inhomogeneous on small scales. In order to use linear 
cosniological perturbation theory one must therefore filter out the small scales, by smearing each 
perturbation over a region of size > 100 Mpc. The same is true at epochs in the relatively recent 
past, except that the comoving filtering scale goes down. Only in the truly early universe is the 
universe (presumably) homogeneous on all relevant scales. 

The filtering is done by means of a 'window function' W{Rf, r), which is equal to 1 at r = and 
which falls off rapidly beyond some radius Rf Taking for definiteness the density contrast, the 
filtered quantity is 

(166) 

(167) 
(168) 
(169) 



and its spectrum is 



where 



and 



The filtered dispersion is 



<5(i?/,x)= / iy(i?/,|x'-x|)5(x')d^a;' 

Vs{Rf,k) 
W{Rf,k) 

Vf = J W{Rf,r)d^x 
a'{Rf)= / W{Rf,k)/Vf Vs{k)- 



W{Rf,k)/Vf\ Vs{k) 
e-'^''W{Rf,r)d^x 



(170) 



The quantity Vf is the volume 'enclosed' by the filter. It is convenient to define the associated 
mass M = poVf^ where po is the present mass density. One normally uses M instead of i?y to specify 
the scale, writing S{M,x.) and a{M). 

The two popular choices are the Gaussian filter 
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and the top hat filter which smears uniformly over a sphere of radius Rf 
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r-Rf) 

3 



Vf = 4:TTRj/3 



W{Rf,k)/Vf = 3 



sin(fci?/) cos(A:i?/)\ 
^kRjr (kRj)^) 
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The Gaussian filter is the most convenient for theoretical calculations, but the top hat filter is widely 
used to as a means of presenting data. 

It is useful to write Eq. (170) in terms of the spectrum Sj^ of the primeval curvature. Using a 
Gaussian filter it becomes 
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On the scales fc ^ > 1 Mpc of most interest, the factor in square brackets increases quite strongly 
with fc, so that the entire integrand peaks on the scale k^^ ~ Rf and 



a^{Rf)^Vs{k)=THk) 



k 



(180) 



with k ^ — Rf. Thus, cr{Rf) falls off as the smearing scale Rf increases. 



Figure 4: The dispersion (To(A/), of the hnearly evolved filtered density contrast evaluated at the 
present epoch, is shown for pure CDM and MDM, with COBE normalization. This figure is taken 
from p], which documents the precise input used. 



One can filte r th e other perturbations as well. The filtered peculiar velocity is called the bulk 
flow. From Eq. ( 115 ), the dispersion of a given component of the bulk flow is 
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It too falls off with increasing scale, but not as fast as the density contrast. 



The bottom-up picture of structure formation 

In Figure 4 is shown the prediction for the linearly evolved a{Rf), evaluated at the present epoch, 
for both pure CDM and for MDM. From Eq. (p7|), the pure CDM prediction scales at earlier epochs 
like (l + z)"!. 

The striking thing about the shape of cr{Rf) is that it decreases monotonically as Rf increases. 
This feature implies a 'bottom-up' picture, with structure forming on successively larger scales. The 
epoch Ziii(M), when a significant fraction of the mass in the universe collapses into objects with 
mass M, is roughly the epoch when a{M) = 1. The linear evolution a oa a = {1 + z)^^ given by 
Eq. ( p7| ) then ceases to be valid, but if it had remained valid the present value of a would have been 

(1 + Znl). Thus 

1 + z„i(M) = ao{M) (182) 

where the denotes the present value of the linearly evolved quantity. 

Figure 4 leads therefore to the following general picture. First, a large fraction of the matter in 
the universe collapses into objects with the mass of a small galaxy, or even smaller. Then successively 
large objects form, presumably gobbling up some or most of the smaller objects. At the present 
epoch, collapse is only just beginning on the scale M ~ IO^^Mq which corresponds to large clusters. 



5.2 Constraining the spectrum of the density contrast 

A full comparison of theory with observation requires the use of numerical simulations, to follow the 
process of gravitational collapse which occurs on each scale after cosmological perturbation theory 



OBSERVATION 


NORMALISATION 


REOIJIRED iT« 


COBE 


crT(10") = (1.1 ± .2) X 10--^ 


1.15± .15 


bulk flow 


(719 .37± .07 


.85± .15 


cluster abundance 


erg = .46± .16 


.46± .16 


galaxy correlation 


1.21 ±.07 


E = .95 


quasar abundance 


ai > .24 at z = 4 


> .4 



Table 2: Key data and their implied normalisation, taken from |3^. The right hand column gives 
the value of ag needed to fit each observation with pure CDM and a flat spectrum, except in the 
fourth row which gives E defined by Eq. (186) 



breaks down. It turns out, though, that the linear theory can be applied on a wide variety of scales, 
so that one can obtain powerful constraints on the parameters by considering it alone. Working 
from the top down, some of the most important linear constraints are explained below. 



The large scale cmb anisotropy The COBE data explores scales of order the size of the observ- 
able universe, say 10"^ to lO^Mpc. 

The bulk flow Smearing the galaxy peculiar velocities over a sphere of radius tens of Mpc to 
get what is called the bulk flow, one should (just) be in the linear regime. In principle Q one 
can observe the radial component of the bulk flow, construct the corresponding potential by 
integrating Eq. (|112| ) radially, reconstruct v and finally determine the density perturbation 
(5p(x) from Eq. ( |11C| ). In practice one has to use in addition the number density contrast Sn/n 
of infrared galaxies (for which the most complete data are available), assume that it is equal 



to bjSp/p with a position independent bias factor 6/, and estimate tJj from Eq. (113) 



Galaxy cluster number density The average number density n{> M) of clusters with mass 
bigger than M ^ lO^^M© gives information on a scale of order lOh^^ Mpc. Within linear 
theory one can estimate n{> M) by assuming that the matter in regions of space where 6{M, x) 
exceeds some critical value 5c of order 1 is bound into objects with mass > M. The fraction 
of space occupied by such regions is 



From this assumption Press and Schechter derived the formula Q 

where v = Sc/<j{M) is the number of standard deviations that 5c represents, and 

(fc2(A/)) = a-^{M) r ee-^"^)Vs{k)'^ (185) 



(This formula includes a more or less unmotivated factor 2). An alternative prescription is 
to identify each peak in i5(M, x) higher than 5c with an object of mass > M , which gives a 
roughly similar result. Yet another method, which in principle is superior, is to run a numerical 
simulation, which again gives roughly similar results. 



The shape of the galaxy correlation function The galaxy correlation function Eq. (121) can be 
used to probe the shape of cr(Af) on scales between those explored by the last two items, if 
the bias factor is taken to be scale independent. The result is usefully quantified by giving the 
'excess power' 

E = 3.4-^;^^fy) (186) 
cr(8/?-iMpc ^ ^ ^ 

where the prefactor is chosen so that E — 1 with pure CDM and n = 1. The observed value 
is somewhat higher. 



Figure 5: The ratio of the MDM and pure CDM predictions shown in Figure 4 is shown, as a 
function of the filtering scale r. Also shown is the result of replacing MDM by pure CDM with a 
tilted spectrum. The figure is taken from pj , where the input is fully documented. 



• Quasar number density Given some astrophysics assumptions, the observed quasar abundance 
can provide a lower limit on the density contrast at high redshift. For instance, one estimate 
pot suggests that at z — A, the fraction of mass bound into objects with M > IO^^Mq is at 
least 1 X lO"'^. Using the Press-Schechter estimate with Sc = 1.33 0, |33|, this fraction is equal 



to 2 times Eq. (183), which gives a{M) > .24 at 2; = 4. 



A subjective view of the normalisation of cr{Rf) on various scales by these observations is given 
in the central column of Table 2, with the notation ar = a{rh~^ Mpc). All values refer to the present 
except the last one. 

First let us see how pure CDM with a flat spectrum fares. The value of o-g required to fits each 
observation is given in the right hand column. The striking thing is that all of the observations 
are consistent with this model to within a factor of 2. On the other hand the observations are not 
actually fitted; normalising to COBE the data fall below the prediction quite sharply as one moves 
down through the tens of Mpc regime. 

In Figure 5 is plotted the ratio of the MDM to the CDM prediction, and one sees that it has 
a sharp drop of about the desired magnitude. Also shown is the effect of tilting the pure CDM 
spectrum, by making n < 1. This too reduces the small scale power, but the drop is not sharp 
enough to accomodate the galaxy correlation data quantified by E. 

One can ask what constraints on n and the fraction il^, of HDM are allowed by the data, if both 
are allowed to vary. Two investigations have been performed [^8[ ^ . The second analysis concludes 
that .20 < f^iy < -35, and that n > .70 (no gravitational waves) or n > .85 (gravitational waves as 
in Eq. ( |165| )), and the first suggests even tighter constraints. Similar bounds on n also hold if one 
tries to rescue pure CDM by throwing away the troublesome galaxy correlation data |33, 2| (justified 
perhaps on the ground that it assumes a scale independent bias parameter). 

These bounds on n are extremely significant in the context of inflation as we discuss in Section 

6. 



5.3 Alternatives to the MDM model 

We have seen that pure CDM roughly accounts for observation, but that to get good agreement 
one needs a rather sharp reduction in power as one goes down in scale from ~ 30h~^ Mpc to 
^ 10h~^ Mpc. The MDM proposal is one way of giving such a reduction. What about others? 

Two have been tried so far. One is to invoke a cosmological constant with Q\ ~ .6, leaving only 
ilm ~ .4 in cold dark matter. While it goes in the right direction, this fix does not seem to be as 
successful as MDM, and is far less pleasing from a particle physics viewpoint. The other is to leave 
the transfer function alone, and require that the primeval power spectrum already has the sharp 
drop in power. Using two inflationary fields one can put this (or practically any other feature) into 
the spectrum, but such 'designer' models are quite unmotivated from a particle physics viewpoint. 
Alternatives which have not been tried yet are to use some kind of 'warm' dark matter. 

Instead of fixing the CDM model one can contemplate throwing it away. Such a radical step 
seems unlikely to be needed, but if one wishes to take it there are two possibilities. One is to invoke 
topological defects to be the seeds of structure formation; the gauge cosmic strings discussed by 
Vachaspati in his lectures might work if the dark matter is hot, and as with the CDM model there is 
continuity as one goes up in scale from galaxies to the scale ~ Hq^ observed by COBE. The other 
is to use a primeval adiabatic or isocurvature density spectrum, but to give it small scale power so 
that the first objects to form have a mass of, say IMq. This would lead to a radically different 
scenario, with no a priori connection between structure formation and the cmb anisotropy. 



6 Inflation 

It is widely supposed that the very early universe experienced an era of inflation il|, 0, |. By 
'inflation' one means that the scale factor has positive acceleration, a > 0, corresponding to repulsive 
gravity. During inflation aH = d is increasing, so that comoving scales are leaving the horizon 
(Hubble distance) rather than entering it, and it is supposed that at the beginning of inflation the 
observable universe was well within the horizon. 

Within the context of Einstein gravity, inflation requires negative pressure p < p/3 (Eq. (|45|)). 
This is achieved if p and p are dominated by a homogeneous, slowly varying 'inflaton field', (pit), 
because one knows from standard field theory that 

p = y+\4>^ (187) 

p = -F + i^^ (188) 

The condition for inflation is clearly 0^ < V . An alternative is to modify Einstein gravity, with or 
without invoking a scalar field, but in that case one can usually redefine the metric so as to recover 
the Einstein gravity during inflation. Until the end of this section I focus on the Einstein gravity 
case. 

It is natural to assume that the field is sufficiently slowly varying that the condition p < p/S 
is amply satisfied, so that the potential dominates p. One then has p ~ —p, which from Eq. ( ^ ) 
makes p practically constant, 

\H-^p/p\ = l+p/p^l (189) 
6.1 General features of inflation 

The most interesting prediction of inflation concerns the adiabatic density perturbation and the 
gravitational waves, which at some level it inevitably generates through vacuum fluctuations. Before 
examining them let us look at the general picture. 

Chaotic initial conditions? 

First let us ask how inflation is supposed to begin, and more generally what happened before it. A 
simple proposal, and the only one that has been explored in detail, is Linde's proposal of 'chaotic' 
initial conditions At some initial epoch a patch of the universe which includes our own as a 
small part, is roughly homogeneous and expanding roughly isotropically. Presumably this initial 
epoch occurs at the Planck scale, with energy density and expansion rate p ^ H ^ ™pii ^^^id it is 
reasonable to suppose that the size of the patch is of order H~^. 



Inflation is supposed to be achieved in the following way. The energy density comes from the 
various scalar fields existing in nature. The initial values of these fields are random, subject to the 
constraint that their potential V{(l), tp, ..) is of order rripi (it is supposed that their other contribution 
to p, coming from the spacetime dependence of the fields, is of the same order rather than much 
bigger.) The inflaton field (j) is distinguished from the non-inflaton fields "ip, ■■■ by the fact that the 
potential is relatively flat in its direction. Before has had time to change much, the non-inflaton 
fields quickly adjust themselves to minimise the potential at fixed 0, after which inflation occurs as 
(j) rolls slowly down the potential. 

An alternative might be to suppose that the universe is initially everywhere roughly homogeneous 
and isotropic, with positive curvature so that there is no boundary. The chaotic hypothesis is more 
flexible in regard to model building, but otherwise the outcome is much the same. 

The = 1 prediction 

We now come to the first virtue of inflation. From Eq. (^), the time dependence of the density 
parameter is given by 



Without inflation, the patch that we are considering will typically either collapse in a Hubble time 
or so (f2 oo) or will become practically empty (fl — > 0). With inflation, in contrast, is driven 
towards 1. Since fl has the same value when observable universe leaves the horizon as it does at 
present, this leads to the expectation that the present value of U, is very close to 1. This conclusion 
might be avoided if the observable universe leaves the horizon near the beginning of inflation, which 
as we remark in a moment is difficult to arrange, or if the initial value of f2 is rather finely tuned 
to zero Whether such escape routes are actually viable is not yet clear, but rightly or wrongly 
the 'prediction' f^o = 1 is widely regarded as one of the firmest consequences of inflation. 

To summarise, we need inflation so that the observable universe is sufficiently old, and is yet 
sufficiently dense to form structure. Naturally, the thought suggests itself that non-inflating parts 
of the universe are common, but simply too inhospitable for us to inhabit. 

Let us note, f or fu ture reference, that after has been driven to a value fairly close to 1, the 
expectation Eq. ( |189|) that p is slowly varying on the Hubble timescale implies the same for H, 
leading to an approximately exponential growth of the scale factor, a cx e'^*. We shall see that to 
have a density perturbation of the right magnitude, p at the end of inflation must be less than about 
(10^^ GeV)'*, so if exponential inflation begins near the Planck scale it must last a very large number 
of Hubble times. This in turn would imply that the epoch at which the observable universe leaves 
the horizon occurs long after the beginning of inflation, because as we shall see this epoch typically 
occurs only 60 or so Hubble times before the end of inflation. 

Perfect classical homogeneity 

So far we need only be dealing only with a patch that is roughly homogeneous and isotropic. As in 
Section 1, the extreme homogeneity and isotropy of the observable universe is 'explained' by asserting 
that it is a very small part of the whole, and that the whole has a non-fractal nature.^ However, 
inflation makes this explanation particularly attractive for two reasons. First, the statement that 
the observable universe is a very small part of the whole patch follows from the exponential growth 
of a combined with the constancy of H, provided that the observable universe leaves the horizon at 
least a few Hubble times after inflation begins 

Second and most crucially, we are free to go all the way and assert that the observable universe 
is absolutely homogeneous at the classical level, or in other words that every Fourier mode of every 
field is in the vacuum state. In fact, not only are we free to make this statement, it is mandatory 
on each scale well before horizon exit. The reason is that even one particle per momentum state on 
such scales would spoil inflation, by dominating the energy density j2j . 

^'^The inflationary 'explanation' of homogeneity and isotropy is often stated incorrectly. Namely, it is said that 
inflation 'smooths out' or even 'stretches' the observable universe. This seems to me to be very misleading. Before 
a scale leaves the horizon during inflation the expansion of the universe is not important, and we can think of 
field oscillations as particles. On the reasonable assumption that interactions are unimportant the the occupation 
number of each momentum state is time independent. After the scale leaves the horizon the field inhomogeneity on 
that scale is frozen in because causal processes cease to operate. This means that (except perhaps for rotational 
velocity perturbations, which would decay anyway) the observable universe is no more and no less homogeneous at 
the beginning of inflation than it is now. 




(190) 



If there is no inhomogeneity at the classical level, how is it generated? The answer, as we discuss 
in the next section, is that an inhomogeneity of well defined magnitude is generated in a perfectly 
natural way by the quantum fluctuation of the inflaton field. 



Avoiding unwanted relics 

Without inflation, thermal equilibrium starts at T ~ 10^^ GeV. Practically every particle is abun- 
dantly produced, and there can be many phase transitions creating defects. Some of these relics 
are definitely unwanted because they would be too abundant, notably examples being magnetic 
monopoles produced at T > 10^^ GeV and relic thermal gravitinos. Such relics can be avoided by 
inflation, because the temperature after inflation (reheat temperature) can be made low enough that 
the offending objects are not produced. Of course, wanted relics might also be killed; in particular 
structure forming cosmic strings need to have energy per unit length of ^ 2 x 10^® GeV, and are 
ruled out by most inflationary models (but see [l^). 



There are two simple reasons for believing that the early universe started with an era of inflation. 
First, without inflation a given patch either collapses (fl oo) or becomes practically empty 
(fi 0) within a Hubble time, unless its density parameter is very finely tuned to the value 17 = 1. 
In contrast, inflation drives 17 towards 1 starting with an arbitrary initial value. Second, inflation 
can postpone the epoch of thermalisation of the universe, so that the temperature is too low to 
produce cosmologically disasterous relics. 

In addition, there is the more complicated issue of homogeneity and isotropy. In order for 
inflation to occur, the energy density of the inflating patch must not be dominated by small scale 
inhomogeneities. In fact, on scales much less than the horizon (Hubble distance) even one particle 
per momentum state would be enough to spoil inflation; in other words the inflationary hypothesis 
requires that the universe is perfectly homogeneous at the classical level, on very small scales. As 
I now explain, the quantum fluctuation on a given comoving scale then generates a well defined 
inhomogeneity and anisotropy, which can be regarded as classical once the scale leaves the horizon, 
and can explain the perturbation the inhomogeneity and anisotropy of the observable universe. 

6.2 The spectrum of the density perturbation 

Inflation predicts a well defined adiabatic density perturbation in the following way. Well before 
horizon entry each Fourier mode of the inflaton field is described by quantum free field theory in 
(practically) flat spacetime. Its mean square quantum fluctuation can be calculated in a straight- 
forward way, being analogous to the Casimir effect for the electromagnetic fleld which has been 
verified in the laboratory. Corresponding to the fiuctuation is a well defined adiabatic density per- 
turbation. Primordial gravitational waves are also generated, through the quantum fluctuation of 
the gravitational wave amplitude. In this section I explain concisely how to do these calculations. 

The slow roll conditions 

A homogeneous scalar field 0(x, t) with minimal coupling to gravity has the equation of motion 



Practically all of the usually considered models of inflation satisfy three conditions, usually 
termed the 'slow roll' conditions. Although the calculations can done without the slow roll conditions, 
they become much more complicated. 

The first slow roll condition is that the motion of the field is overdamped, so that the 'force' V' 



Summary 



(l} + 3H<j) + ¥'{(/)) = 



(191) 



balances the 'friction term' 3iJ0, 




(192) 



The second condition is that 




(193) 



which means that the inflationary requirement (p^ < V is well satisfied and (assuming f7 = 1) 



^ (194) 



These two conditions imply that H is slowly varying, and that the scale factor increases more or 
less exponentially, 

a (X e"* (195) 

The third condition is that 

\v\ < 1 (196) 

where 

ml, V" , , 



It can be 'derived' from the other two by differentiating the approximation Eq. (|5^) for and noting 
that consistency with the exact expression Eq. ( ^Sf ) requires (?i ^ is satisfied. However there is 
no logical necessity for the derivative of an approximation to be itself a valid approximation, so this 
third condition is logically independent of the others. 

It should be noted that the first slow-roll condition is on a quite different footing from the other 
two, being a statement about the solution of the field equation as opposed to a statement about the 
potential that defines this equation. What we are saying is that in the usually considered models 
of inflation, one can show that the first condition is an attractor solution, in a regime typically 
characterised by the other two conditions, and that moreover reasonable initial conditions on (j) will 
ensure that this solution is achieved well before the observable universe leaves the horizon. 



Classical equations 

To define the perturbation Scj) of the inflaton field, one has to choose a slicing of space-time into 
spacelike hypersurfaces, just as was the case for the perturbations Sp, Sp and SH that we studied 
earlier. Since the inflaton field is supposed to dominate the energy-momentum tensor, the momentum 
density vanishes if its spatial gradients vanish In other words, Scl) vanishes if the hypersurfaces 
are chosen to be comoving! 

In this situation one can proceed heuristically in the following way ^. First one notes, from 
Eq. (p^, that in the extreme slow roll limit p/p —1 (corresponding to p ^ 0), the curvature 
perturbation TZ of comoving hypersurfaces becomes infinite, if the comoving density perturbation 
is finite. One should therefore define the inflaton fleld perturbation 6(j) on a family of hypersur- 
faces which remains undistorted in the slow roll limit. The time displacement 6t of the comoving 
hypersurfaces from the undistorted ones is determined by the fact that Scj) vanishes on the latter. 



(198) 



Finally, one can show either by brute force [ pT| or by an elegant geometrical argument j^] that 
the distortion of the comoving hypersurfaces caused by this time displacement implies a curvature 
perturbation TZ = —H5t, which in terms of the field perturbation is 



7^ = H6(j>/<p 



(199) 



Alternatively one can derive this equation without any mention of the spatial geometry, working 
entirely with the density perturbation plus the definition Eqs. (|9l|) and (BSh of TZ p9[. 



Perturbing the inflaton fleld equation Eq. (191) leads to 



+ V'' 



= 



(200) 



Until a few Hubble times after horizon exit Eq. ( p^ ) ensures that V" can be dropped, so that the 
equation becomes 



(^0k)"+3ff(<50k)' + 







(201) 



A rigorous treatment H, M is to define Scfi on hypersurfaces with zero curvature perturbation, 



Eq. (199) then being an exact expression for the curvature perturbation of the comoving hypersur- 
faces |pl"|. One can show that to first order in the cosmological perturbations the field equation for 



is 

u2 



< + (fc^ - z"Mk) = (202) 
where u — aS(j), z — a(j)/H and a prime denotes a{d/dt). Finally, one can show that this equation 



reduces to Eq. (201) if the slow roll conditions are satisfied. (The correction can be calculated 



explicitly, and has been shown to be small [[16|.) 
The quantum fluctuation 

Well before horizon entry 50k is a massless field living in (practically) flat space-time, since its 
wavenumber k/a is much bigger than the Hubble parameter H. It can be quantised in the usual way 
so that its quantum state is labelled by the number of inflaton particles present with momentum k, 
and we are assuming that there are no particles, corresponding to the vacuum state. 
Working in the Heisenberg representation, (/)k is associated with an operator 

Mt) = Mt)ak + wUt)alj^ (203) 
The annihilation operator Sk satisfies the commutation relation 

[aki,4 J ='5ki,k2 (204) 
and the field satisfies the commutation relation, 

d - 

[0(xi,i),— (?i(x2,i)] = iS^iaKi-ax2) (205) 



m^(5^(xi-X2) (206) 



As a result the functions Wk are given by 

Wkit) = a-3/2(2fc/a)-i/2g-^(x+(fct/a)) ^207) 

The phase factor x is arbitrary, and it has negligible variation on the timescale a/k. The vacuum 
state is the one annihilated by dk, so the vacuum expectation value of the field perturbation is 

(|^0kP) = \wk\' (208) 

To extend these results to the epoch of horizon exit and beyond, we have to accept the validity 
of free field theory in curved space-time. All we need to assume is that there is a Heisenberg picture, 
in which operators satisfy the classical equations of motion and state vectors are time independent. 
Then (^k continues to satisfy the field equation Eq. ( p00| ), and E qs. (p03| ), ( |204| ) and (|20|) stiU hold. 



where Wk is the solution of the field equation reducing to Eq. (207) well before horizon entry. As 
one easily checks, the required solution is 

A few Hubble times after horizon exit, the vacuum expectation value is therefore 

(l^0k|') = ^ (210) 
A measurement of the 0k 's will yield random phases, and a distribution of moduli whose disper- 



sion is given by Eq. (21C). Accordingly the spectrum of the inflaton field, defined by Eq. 



given a few Hubble times after horizon exit by 

-P'J'ik) = ^ (211) 
^ Ztt 

Since H is slowly varying on the Hubble timescale it can be written, a few Hubble times after horizon 
exit, 

r'J'ik) = ^ (212) 

^ ZTT 



where the star denotes the epoch of horizon exit, and is evaluated a few Hubble times after 
horizon exit. 

The spectrum of TZ is given by 

VU' - "^VT (213) 

In Section 2 we learned that 7?,k is constant after horizon exit, so Vn remains constant even though 
H/(j) and might vary separately. It follows that as long as the scale is far outside the horizon, 

^ — ^ (214) 
27r0» 

where the star denotes the epoch of horizon exit. Using Eq. (|57|), this leads to ||2^ 



5l{k) = (215) 
75 m 



PI 



Primordial gravitational waves? 



The gravitational wave amplitude hij is dynamically equivalent to a minimally coupled, massless 
scalar field ?A — {mjn/lGTrY^^hij . Each non-zero Fourier component has the same vacuum fluctua- 
tion as a massless field, with the above conversion factor. Thus the spectrum of the gravitational 



waves is given by Eq. (211) 



(216) 



Putting this expression into Eq. (15C) and dividing it by Eq. (|146|) gives the ratio Eq. (|152|) of the 



gravitational wave and density contributions to the cmb anisotropy. 



The classicality of the perturbation 

The quantum fluctuation of a field in flat spacetime cannot be regarded as classical. The reason is 
that each Fourier mode is equivalent to a harmonic oscillator in its ground state; the position and 
momentum of such an oscillator cannot be accurately measured simultaneously (unless the measured 
results are far out on the tail of the probability distribution), and as a result the position cannot be 
well defined over an extended period of time. Returning to the case of the field, this means that a 
Fourier component cannot be defined with high accuracy over an extended period of time. 

After horizon exit though, a given Fourier component is no longer equivalent to a harmonic 
oscillator and one can show that it can have a well defined value In this sense it is classical, 

though of course one still has the 'Schrodinger's cat' problem of how and when it is supposed to 
have acquired this well defined value, drawn from the quantum probablity distribution. 



The inflationary energy scale 

Assuming that the adiabatic density perturbation dominates the cmb anisotropy, we saw earlier that 



gives Sh — 2.0 x 10 on the Hubble scale. Eq. (211) therefore gives the energy scale at the epoch 



when the observable universe leaves the horizon during inflation, 

V,'/^ = 7.3 X IQie GeV x ( \ ^'^ (217) 

^ V2.0xl0-V 

Since ei <C 1 we conclude that V^l^ is less than a few times 10^® GeV. Since the energy density 
decreases with time, this is also an upper limit on the energy density at the end of inflation, and on 
the reheat temperature. 



6.3 Entering and leaving inflation 



In order to work out the observational consequences of Eq. (217) for particular inflationary models, 
we need to do a little more work. 

First, we need to ask how inflation ends. In the chaotic picture, two possibilities exist concerning 
the minimum into which the non-inflaton fields fall. The simplest possibility is that it corresponds to 



the true vacuum; that is, that the non-inflaton fields have the same values as in the present universe. 
Taking the potential to be quadratic for simplicity, it then has the form 



V = ^m^<j)^ (218) 

Inflation ends when the inflaton field starts to execute decaying oscillations around its own vacuum 
value (zero in the above example). This typically occurs when the slow roll conditions fail, 

max{eo„d, |?7ond|} ^ 1 (219) 

The hot Big Bang ensues when the vacuum value has been achieved and the decay products have 
thermalised. Models of this type have been widely explored 

The other possibility |jl|, |, |l^, [l8[ |l|| is that the minimum at large fixed <f> corresponds to a false 



vacuum, which dominates the 0-dcpendent part of the potential. In the quadratic case the potential 
is then 

V ^Vq + ^771^0^ ~ Vb (220) 

where Vq is the false vacuum energy density. In this case inflation ends when the false vacuum is 
destabilized, as (j) falls through some critical value (pc- 

Of course one can envisage more complicated possibilities. One is that an early epoch of inflation 
gives way to a hot universe, which finds itself in a false vacuum which is stabilized by the finite 
temperature. The vacuum energy density Vq dominates when the temperature falls below Vq^^ , and 
infiation occurs for a more or less brief era as the field rolls slowly away from its minimum, ending 
when it starts to oscillate around its true vacuum value. This way of beginning inflation is called 
'new inflation' [Q. Apart from being more complicated, it is also even more difficult to implement 
in the context of sensible particle physics than the other two proposals. 



Reheating 

In any model, the universe (or at least the constituents dominating the energy density) must achieve 
thermal equilibrium at some point after inflation. Because the first models of infiation invoked a 
thermal entry, this event is usually called 'reheating'. The corresponding 'reheat temperature' is the 
biggest temperature ever achieved after inflation and plays a crucial role in cosmology. Unfortunately, 
no reliable estimate of it is known at present. If inflation ends with a phase transition, reheating is 
expected to be prompt. If it ends with the inflaton flcld oscillating, reheating may be long delayed 
because the couplings of the inflaton field to other fields is then typically quite weak. 



The epoch of horizon exit 

A first requirement for an inflationary model is that the observable universe should be within the 
horizon at the beginning of inflation. We therefore need to know the epoch when a given scale leaves 
the horizon, something which is necessary also to calculate the inflationary perturbations. Denoting 
it by a star, this epoch is given by 

a^H^ = k (221) 

The epoch of horizon exit is therefore related to the present magnitude of the scale in Hubble units 
by 

aoHo^aoHo ^222) 
k CL^H^ 

Let us denote the end of inflation by a subscript 'end' and the epoch of reheating by 'reh', assuming 
matter domination during the era (if any) between these two epochs. Let us also assume that after 
'reh' there is radiation domination, until the epoch 'eq' at which the dark matter density becomes 
equal to that of the radiation. Throughout the history of the universe the radiation energy density 
is proportional to a~*, that of the matter is proportional , and the total is proportional to . 
It follows that 







(223) 
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e'^' f^) ei) (224) 
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where por = (acq/ao)po is the present radiation energy density and A^* is the number of Hubble 
times between horizon exit and the end of inflation, 
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It foUows that 
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This equation relates the three energy scales V^^^, KmdL^'^'-L^tsh ■ ^^'^^ ^^"^ scales are related 



by another expression for A^,, which follows from Eqs. (226), (192), (193) and (194) 
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The biggest scale that can be explored is roughly the present Hubble distance, ao/k 
3000/i~i Mpc. I shall refer to the epoch of horizon exit for this scale as 'the epoch when the observable 



Ho' 



universe leaves the horizon', and denote it by the subscript 1. As displayed in Eq. (217), the COBE 



observations require that Vi^'^ < IQi^ GeV, with the equality holding in most models of inflation. 



Also, Eq. (228) gives in most models of inflation V^^^ ~ Knd'- reheating is prompt, we learn 



.1/4 



that in most models the observable universe leaves the horizon about 62 e-folds before the end of 
inflation. If reheating is long delayed A^i could be considerably reduced, being equal to 32 for the 
most extreme possibility of V^^^^^ ^ 1000 GeV (corresponding to reheating just before the electroweak 
transition which is presumably the last opportunity for baryogenesis). For most purposes, however, 
one needs only the order of magnitude of A^i. 

The smallest scale on which the primeval perturbation can be probed at present is around 1 Mpc, 
and one sees that this scale leaves the horizon about 9 Hubble times after the observable universe. 
Discounting the case of very late reheating, we conclude that in the usually considered inflationary 
models, scales of cosmological interest leave the horizon 50 to 60 Hubble times before the end of 
inflation. 



The spectral index 

The spectra l ind ex n o f the d e nsity pertu rbat ion is obtained by differentiating Eq. ( |215| ) with the 
aid of Eqs. (|l9|), (^, (|l9|), (|2l|) and (|95|) 



n = 1 + 2771 - 6ei (229) 
For definiteness I have taken to be the epoch when the observable universe leaves the horizon. 

6.4 Specific models of inflation 

Of the many viable inflationary potentials that one could dream up, only a few are reasonable from 
a particle physics viewpoint. Some of them are described now. 



True vacuum chaotic inflation 



If the potential is V (x 0", Eq. (21£) shows that at the end of inflation (pcnd — a(mp;/87r)i/'^. From 
Eq. (|22|) with A^i = 60, 0? ~ 120aTO2,;/87r. This leads to 1-n = (2 + a)/120 and R = .05a, leading 
to the following relation between the gravitational wave contribution and the spectral index 



i? = 6(1 - n) - 0.1 



(230) 



For a = 2, 4, 6 and 10 one has l~n— .033, .05, .067 and .1, and R — .10, .20, .30 and .50. Lazarides 
and Shaft [Q have shown that a potential of this kind follows from a class of superstring-inspired 
gauge theories, with an index a varying between 6 and 10. 



False vacuum chaotic inflation 



The simplest model of false vacumn inflation is [|T6[ |[ |l^ |T^, |l9| 

V{^, <f) - iA(^2 - M^ f + im202 + Ia W (231) 

The couplings A and A' are supposed to be somewhat less than unity, and in the numerical examples 
I set A = A' = 0.1. For cj)^ > = AM^/A', the potential for the ■0 field has a local minimum at 
"0 = 0, in which the field is assumed to sit (except for the quantum fluctuation). Inflation can occur 
as <j) slowly rolls down its potential 

V{(j)) = ^AM^^ + ^m202 (232) 

We will suppose that while observable scales leave the horizon the first term dominates, since in the 
opposite case we recover the 0^ potential already considered. When falls below (t>cnd this local 
minimum becomes a local maximum, and a second order phase transition to the true vacuum occurs 
which ends inflation and causes prompt reheating. 

The requirement that the first term dominates means that 

1^ « > (-) 

Of the parameters e and -q which are required to be small, the second is independent of 0, 

V = (234) 



where 



2 _ mil m2 



The ratio X must therefore be significantly less than 1. From Eq. (|22§|), 



01 = VAAV^^e^i" (236) 
Consistency with Eq. ( ^33| ) requires roughly 77 < .1. The other small quantity ei is given by 

1 A Stt 



ei 



2 A' m% 



^j2^2g2A'i„ ^237) 



and Eq. (233) requires ei <C rj. It therefore follows that the spectral index n is bigger than 1 in the 
two-scale model. 

Setting A = A' = .1 and imposing the COBE normalisation 6h = 1.7 x 10~^ determines all of the 
parameters in terms of m. For m — 100 GeV one has M = 4 x 10^^ GeV leading to 77 = 10""* and 
Ci = 10^^^. The gravitational waves are absolutely negligible, and the spectral index is extremely 
close to 1. The maximum value of m permitted by Eq. ( ^33| ) is roughly m = 10^'^ GeV, giving 
M = 2 X 10^^ GeV, 77 — .07 and ei = lO^'^. The gravitational waves are still negligible, but n — 1.14, 
significantly bigger than 1 . 

So far we have considered a single scalar field ip, possessing the discrete symmetry ip — ?/). 
When inflation ends, domain walls will form along the surfaces in space where ip is exactly zero, so 
to make the model cosmologically viable one would have to get rid of the walls by slightly breaking 
the symmetry. However, it is clear that one can easily consider more fields, with continuous global or 
gauge symmetries. In particular, if ip is complex one can use the same potential with the replacement 
■0^ \ip\'^. This leads to exactly the same inflationary model, but now global strings are formed 
instead of domain walls. 

The case m ^ 100 GeV and M ~ 10^^ GeV is particularly interesting because these are rather 
natural scales in particle physics. One possibility is to identify ip with the Peccei-Quinn field, and 
it might also be possible to realise the model in the context of superstrings |lq] 



Natural Inflation. 



The 'natural inflation' model Q has the potential 

1/(0) = A4(1±cos(0//)) (238) 
where A and / are mass scales. In this model 

e = r-^tan^^ (239) 
V = -r^^ f 1 - tan^ ^ ) (240) 



2/ 

where r = (167r/^)/TOp;. Since e — rj ~ r^^ one must have r 3> 1 to be in the slow-roll regime 
anywhere on this potential. This requires that / is at the Planck scale, suggesting a connection with 
superstrings which has been explored in jSOt . 

In this model inflation has to begin with close to zero. If there is a preceding thermal era, 
one might expect (j) to be chosen randomly so that a value close to zero is certainly allowed. With 
or without a thermal era, additional flelds have to be involved to have a chaotic beginning at the 
Planck scale. 

Inflation ends with an oscillatory exit at the epoch given by Eq. ( pl9D , tan(0ond/2/) — r^^^ 
which is of order 1. Using this result, Eq. ( |228D with Ni = 60 gives (/)i/2/ ~ exp(— 60/r), leading to 

ei = (241) 
r 

771 ~ -l/r (242) 

l-n ~ 2/r (243) 

Thus, natural inflation makes the gravitational waves negligible but tends to give a significantly 
tilted spectrum. The observational bound n > .7 implies that r > 6. 

i?^ inflation 

The Lagrangian density for Einstein gravity is just — (mp;/167r)i?, where R is the spacetime curvature 
scalar. Perhaps the simplest model of inflation is to replace Rhy R + R^ / (6M^). By redefining the 
metric one can then recover Einstein gravity, at the expense of modifying the rest of the Lagrangian 
and adding a completely new field (j) with a potential 



327r 



1 — exp 



1« \ 1/2 



3m|,; 



(244) 



In the regime >^mpi, the potential V{4>) satisfies the slow- roll conditions so that infiation occurs. 
The non-gravitational sector is irrelevant during inflation, and if the R^ term quickly becomes 
negligible afterwards we have a model in which infiation is driven by a scalar field, without any 
modification of gravity. What has been achieved is to motivate the otherwise bizarre form of the 
infiationary potential. 

In the regime > mpi where they are small, the parameters appearing in the slow-roll conditions 

are 



(245) 
(246) 

There is an oscillatory exit to infiation at the epoch given by Eq. (^19D, 0cnd ("^pi/STr)^/^. From 




Eq. ( |228| ) with TVi = 60, 

^1 - 5^ (247) 

leading to rji ~ —.02 and ei ~ lO^**. Thus the spectral index is n = 0.96, but the gravitational wave 
contribution is completely negligible. 

Justification for this model of infiation has recently been claimed in the context of superstrings 

ill 



Extended inflation 



'Extended inflation' models are based on modifications of Einstein gravity of the Brans-Dicke type 
p2| . After recovering Einstein gravity by redefining the metric, one typicaUy has an infiationary 
scalar field, with the potential 

V = Voexp{ J^^<j)\ (248) 




Irrespective of when inflation ends, one finds —Ug = 1 — n(— 2/p), and the 'canonical' gravitational 
wave ratio R ~ 6(1 — n) that we discussed earlier. 

Infiation in these models ends with a first order phase transition (caused, in the Einstein gravity 
picture, by the modification of the non-gravitational Lagrangian induced by the metric transfor- 
mation). The bubbles that it generates must not cause an unacceptable microwave background 
anisotropy, which generally leads to the constraint n < 0.75 |3^, in confiict with the observa- 
tional bound n > .85 that we noted in Section 5. One needs to make these models quite complex in 
order to avoid this difficulty jHJ, |l^ . 

6.5 Summary and prospects 

Particle physics, especially in the context of superstrings, is beginning to hint at several viable models 
of inflation. Each of them makes a distinctive prediction regarding the tilt of the density perturbation 
spectrum, which will allow one to discriminate between them in the near future provided that some 
variant of the CDM model continues to fit the observations. One will then, for the first time, have 
an observational window on energy scales not far below the Planck scale. 
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